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Abstract 

Q-ij A manifestly gauge-invariant hamiltonian formulation of classical 

electrodynamics has been shown to be relativistic invariant by the con- 
struction of the adequate generators of the Poincare Lie algebra |Q. 
^ • The original formulation in terms of reduced distribution functions for 

the particles and the fields is applied here to the case of two charges 
interacting through a classical electrodynamical field. On the other 
hand, we have been able in previous work to introduce irreversibility 
at the fundamental level of description Q by reformulating field the- 
ory without self-energy parts by integrating all processes associated 
with self-energy in a kinetic operator, while keeping the equivalence 
c/2 ' with the original description. In this paper, the two approaches are 

combined to provide a formalism that enables the use of methods of 
statistical physics j| to tackle the problem of the divergence of the 
self-mass. Our approach leads to expressions that are finite even for 
point-like charged particles: the limit of a infinite cutoff can be taken 
in an harmless way on self consistent equations. In order to check 
our theory, we recover the power dissipated by radiation in geometries 
. where the usual mass divergence does not play a role. 
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1 Introduction 



The derivation of an equation of motion of an electron that includes its 
reaction to the self-field has been initiated by Abraham and Lorentz hundred 
years ago and is still a controversal matter Q , || . The main problem is the 
presence of divergences associated with point-like charged particles. A way 
of removing them has to be devised without entering in trouble with the 
special theory of relativity (see Ref. || for a recent review and a relevant 
bibliography). The derivation of the self- force based on energy conservation 
|?J avoids that problem: the power emitted in the radiating field is due to 
the work of the radiative reaction force. 

In the usual derivation, a given motion is prescribed to the charge and 
the potentials of Lienard-Wiechert associated with it are computed. The 
self-fields are then derived and their expression used to get the reaction on 
the motion of the charged particle, leading to an (infinite) term interpre- 
tated as a self-mass. But an infinite self-mass prevents the acceleration, 
hence the paradox. An approach that uses a finite expansion in the charge 
does not allow to correct the situation. In a theory without divergence, a 
natural cut-off for the wave numbers, proportionnal to the inverse of the 
classical radius of the electron, should appear but that property cannot be 
obtained in a simple expansion in the charge since the cut-off value should 
be proportionnal to the square of the inverse of the charge. Moreover, the 
use of an effective frequency dependent mass M(u>) @ leads to an equation 
that is not free of runaway solutions. 

In different contexts, methods of statistical mechanics || enable, through 
resummations of formally divergent class of diagrams, to get relevant finite 
physical results. Here we present a formalism adapting those methods for 
problems in classical electrodynamics. 

A manifestly gauge-invariant hamiltonian formulation has been devel- 
oped for systems composed by point particles and fields, the state of which 
is now described by reduced distribution functions The dynamical vari- 
ables are the positions and mechanical momentums for the charged point 
particles and the transverse components of the electric and magnetic fields. 
Potentials are absent in that formulation. The Coulomb interaction takes 
into account the longitudinal part of the electric field. A generalized Liou- 
ville equation for the reduced distribution functions is derived. It provides 
a statistical description that takes into account the Lorentz force between 
the particles and the Maxwell equations for the fields. The formalism for- 
mally ressembles a statistical description of charged particles in Coulomb 
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gauge but with a different interpretation and the guarantee of satisfying the 
principle of special relativity. The relativistic invariance is proved by the 
explicit construction of the generators of the Poincare Lie algebra. 

While Balescu and Poulain have developed their formalism for an arbi- 
trary number of particles, described by reduced distribution functions, we 
can apply it as such in the simplest case of two charged point particles. They 
thus interact through the Coulomb interaction and the classical transverse 
electrodynamical field (electric and magnetic). An alternative possibility is 
to consider a single charge in interaction with a Coulomb potential (due for 
instance to an infinitely massive particle) at the origin of the coordinates but 
the translation invariance is then immediately broken. At the final stage, for 
the sake of interpretation, we will consider the limit of our expressions when 
one of the two particles becomes very massive. Working with two particles 
avoids the consideration of an external force to accelerate the particles: the 
relativistic and gauge invariance is therefore preserved. A prescribed mo- 
tion for the charged particle is therefore avoided. The consideration of an 
incident transverse field is also relevant to the problem but is not treated 
here. 

The Balescu-Poulain formulation seems therefore an adequate starting 
point to deal with classical electrodynamics thanks to its intrinsic properties: 
namely relativistic invariance, explicit gauge invariance. The formulation is 
statistical: the particles and the modes of the field are described by distri- 
bution functions. The distribution functions associated with the particles 
can be spatially extended. A particle does not interact with the electric 
longitudinal field its generates: the Coulomb interaction is considered only 
between different particles. The present paper starts with the results of the 
last section of the paper of R. Balescu and M. Poulain [jy]. 

A theory of subdynamics has been introduced thirty years ago by the 
Brussels group (see e.g. |§, ||) for a dynamics provided by the Liouville-von 
Neuman equation. A setback of that approach is a limitation on the class 
of possible initial conditions since they have to belong to the subdynam- 
ics. To avoid the trap, we have introduced the so-called single subdynamics 
approach || based on the existence of self-energy contributions to the dy- 
namics. In that way, we obtain a reformulation of field theory that excludes 
self-energy contributions in the dynamics. However, being able to accom- 
modate also initial conditions outside the scope of the original dynamics: 
our dynamics is larger that the initial formulation. Since the formal prop- 
erties of the subdynamics do not depend on a particular realisation of the 
operators, we have picked up all the formal properties without a need to 
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redemonstrate them. The so-called single subdynamics approach has been 
illustrated first in a quantal non relativistic framework @, Q, Q, gH, 
|^| . All the claimed formal properties have been subject to an explicit check 
for exactly soluble models. These basis ideas can also be applied in the 
present classical relativistic context. 

The adequate way of dealing with the self-field is provided automati- 
cally by the single subdynamics approach. The dynamics is first extended 
to be able to distinguish the self-field contributions from the other. A sub- 
dynamics, inspired by the formalism developed at Brussels || enables the 
obtention of dynamical equations of motion in which the self-field does no 
longer appear. It has been proven [13] that the description exactly contains 
the original description and that the effects of the self-field are now taken 
into account in the new generator of motion. The relevant subdynamics in- 
corporates all the features of usual CED. That description therefore includes 
not only the original dynamics but could also include a more general class of 
initial conditions, enabling the inclusion of irreversibility at a fundamental 
level. 

Here, we deliberately restrict ourselves to the derivation of the closed 
irreversible evolution equations for the interacting charged particles, in the 
absence of incident field from an outside source, and to the obtention of the 
emitted fields (velocity and acceleration fields) at the lowest non-vanishing 
order. 

Our paper is structured as follows. First, we present the Balescu-Poulain's 
formalism and derive the evolution equations for all reduced distribution 
functions defining the state of the system. The basic idea for constructing 
the single subdynamics in CED is the use of a distinction between real and 
virtual fields (the virtual field forms the self-field). We propose an extension 
of the dynamics suitable for our purpose and the constitutive relations that 
connects the original and extended dynamics are displayed. The elements 
ot the extended dynamics bear a tilde accent. 

The kinetic operator 0, considered in section 3, describes the closed evo- 
lution of the distribution functions that do not involve the self-field. Their 
elements are evaluated from the corresponding vacuum-vacuum elements 
of the subdynamics operator E(t). The first non- vanishing contribution 
appears at the second order in the interaction with the transverse fields, 
without considering, in the first step, the influence of the Coulomb interac- 
tion between the charged particles. The various steps of the derivation are 
illustrated and the final expression for is given in Appendix B. All the 
elements are known to examine the putative second order mass correction 
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for the charged particles, that is found to vanish. 

A non- vanishing contributions to the kinetic operator, reflecting the pres- 
ence of the effect of the transverse self-field, requires to consider either a 
non-vanishing incident transverse field, either a Coulomb interaction be- 
tween the charged particles or either the mutual influence of the transverse 
emitted field: the particle has to be accelerated to receive a radiation reac- 
tion force. 

To get a better insight of the previous result, we take another road in 
§4. The kinetic operator can indeed be also evaluated from the knowledge 
of the self-field determined by the so-called creation operator. The value 
of the self-field at the location of the particle induces its self-interaction. 
Since the equivalence conditions require the equality of the emitted and self- 
field, the creation operator provides us moreover with the expression of the 
emitted field. Correlation-vacuum elements of the resolvent are considered 
for evaluating the elements of the subdynamics. A simple computation 
enables to get explicitly the expression of the common value of the Fourier 
transform of the emitted and self- field. 

To obtain a source of acceleration and to prepare an easy comparison 
with the usal approaches, the first order effect of the Coulomb acceleration 
will also be computed in §5 and §6 from two different ways: the direct con- 
sideration of the kinetic operator and the recourse to the creation operator 
for the self-field. 

The direct computation of the kinetic operator is performed in the next 
section §5 from the vacuum- vacuum elements of the resolvent acting in ab- 
sence of field (field vacuum). All relative orders of the vertices have to be 
considered: the Coulomb interaction can a priori take place before, after or 
between the two interactions with the transverse field. Only the last two cir- 
constances lead to a non-vanishing contribution. Indeed, when the Coulomb 
interaction takes place after the two interactions with the transverse field, 
we receive as factor, as expected, the previous vanishing second order con- 
tribution to the kinetic operator. The computation, although lenghty, is 
straightforward. 

For a consistency check, in §6 we consider the creation operator at first 
order in field-particle interaction and first order in the Coulomb interaction. 
This enables to get the effect of the acceleration, due to the Coulomb in- 
teraction, to the self-field, hence to the retroaction of the emitted field on 
the accelerated particle. From the equivalence conditions, we deduce for 
all points the field emitted during the acceleration of the particle. If we use 
that expression in the kinetic equation, we recover the previous result. From 
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its expression at the localisation of the particle, the power emitted can be 
computed. 

Our expressions are analysed in §7. We consider a situation in which 
the distribution functions of the charged particles are infinitely sharp in 
configuration and momentum space, with an absence of free field. One 
particle is then considered as infinitely heavy and we use the referential in 
which the heavy particle is at rest. In the geometry where the position and 
velocity vectors are orthogonal, the power dissipated by the field due to 
the motion of the light particle can be computed exactly: all integrals can 
be performed. In other geometries, we do not avoid the usual divergence. 
This is natural since our approach contains the usual formalism and no 
resummation has been performed yet. The usual result is explicitly recovered 
as a particular case in small velocities circonstances. Indeed, under the 
equivalence conditions, both theories provide the same equations for the 
motion of the charged particles. 

The finiteness of the theory through resummations is considered in §8. 
We start from the divergent contribution for the self-field at the lowest order 
in the charge. We regulate it by introducing a cut-off function, with a cut-off 
value K c to smoothen the contributions from very high wave numbers. That 
contribution is the first term of a series that can be formally resummed, in a 
self consistent way. We then show that the limit K c —* oo can be performed 
in an harmless way. The effective cut-off resulting from the solutions of the 
non-linear equations is naturally linked to the classical radius of the electron. 
Our final expressions do not admit a simple expansion in the charge, since 
the dependence of the effective cut-off in the charge e is in -j. This result 
is obtained through an analysis that uses several steps. 

Some conclusions and perspectives are considered in the last section §9. 

2 Electrodynamics in the manifestly gauge in- 
variant Balescu-Poulain formalism 

In this section, we define the model for the description of the two charges 
in interaction with the electromagnetic field. We use the approach by R. 
Balescu and M. Poulain [ffl. The only difference is that for the description of 
matter, we do not deal with a reduced formalism but keep the two-particle 
distribution function. Although the transposition is straightforward, we will 
present it in the following subsection, using their notations (and expressions 
whenever possible). 
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2.1 The Balescu-Poulain formalism 



The state of the system is described by a distribution vector J 7 , i.e. by 
a collection of functions describing two different particles and the reduced 
distribution of m field oscillators, describing the transverse field components 
which are the only ones that appear explicitly: 

T = {/ii H (x {1) ,x^;x [l K.-.,X [m] ; k [1] , • • • , kH )} ; m = 0, 1, 2, . . . 

(2.1) 

For m = 0, the system does not contain fields variables. An obvious conven- 
tion in the notation is implicit for m = 0. Here x^> denotes the coordinates 
(<! jP^') °f particle j, and denotes the variables describing a given 
field oscillator associated with the wavevector k^: (j]a , £a , ct = 1, 2) 111 that 
are the action ) and angle variables (£a ) associated with the oscillator 
characterized by the wave number k^l. 

If two mutually orthogonal unit vectors, or "polarization vectors", e°(k) 
associated with a given wavevector k are introduced such that, together 
with the unit vector ¥, they form a right-handed cartesian frame, the elec- 
tromagnetic fields are expressed as follows in these variables. 

E ± (x) = —^ E /d 3 A;fc5e Q (k)7 ? l(k)exp{m[k.x-2< Cf (k)]}, 
(2vr)2 a=h2a= ±i J 

(2.2) 

B(x) = — E /rf 3 fcA;i(-l) Q 'e Q '(k)r/l(k)exp{m[k.x-27re a (k)]}, 
(2vr)2 a=lj2a =±i J 

(2-3) 

where a' = 2 for a = 1 and a' = 1 for a = 2. 

The dynamical functions of the system are described by a set B: 

B = {b n[m] (x« , x® ■ x W , . . . , X H . k [H , . . . , k M )} ]m = 0,1,2,... 

(2.4) 

The average value of an element &n[ m ] of B is calculated by the following 
formula: 

< b n[m] > = f d 3 kW . . . d 3 kH f d \W . . . d \[m] f d 6 x (l) ^(2) 



1 In contrast with |J_5[, the reduction is not performed up to the level of each polarized 
mode, in the same way that reduced distribution functions for the particles are not con- 
sidered to only one component of the velocity. This procedure ensures more easily the 
rotationnal invariance of the treatment. 
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x & 11[m] (^),x( 2 );x [1] ,...,X H ;kW ) ...,kH) 

x fuH , * (2) ; x [1] , • • • , x [m] ; kW , . . . , kH ) . (2.5) 

It is well known that to each generator G of the Poincare Lie algeba 
corresponds an infinite hierarchy of equations describing the transformation 
properties of the reduced distribution functions [10] ] . These equations can 
be written compactly as 

d g T = C G T, (2.6) 

where Cq is a matrix operator. The components of this equation are written 

as 

oo 

9gfu[m] = E < 11M|£ G |11K] > Aim- (2-7) 

m'=0 



The matrix elements entering these equations are obtained as in |T7j] and 
listed below, considering separately the three contributions corresponding 
to the splitting of the liouvillians in three terms, describing respectively free 
particles L^f, free field L^f and interactions L' G . 
For the free particles, we have: 

< llH|£g p |ll[m'] >= 8 mm , (4 (1) +4 {2) ) . (2.8) 

In this work, we consider only the generator corresponding to the time trans- 
lation (g = t, G = H) and get, using Einstein convention for the summation: 

^ = -^4, (2-9) 
oqr 

(i) (i) 
where the velocity vf is connected with the mechanical momentum pf' in 

the usual way (in the units chosen, c = 1 and the div (divergence) of the 

electric field vector is 4-7T the charge density): 

„(i) 

V P = , 2 u) ov - (2 ' 10) 

(rrij +ps ps ) 2 

For the free field, we have: 

m 

< ll[m]|4f|ll[m'] >=6 mm ,J2Lf, (2.11) 

i=l 
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a 



a=l ue ia 



For the interaction, 

< ll[m]|£' G |ll[m'] >= 5,, 



(2.12) 



4. f'^UrS , t' {12) 



i=l 



+ <W, m +i / <?k^ J dj^ (Lt +1] + 4 2[m+11 ) , (2-13) 
where / d^ m+1 ^ stands for 

/ J 7 K+i] . . . = f° 4 m+1] /°° d4 m+1] ^[ m+1] ^ m+1] . . . (2.14) 
J Jo Jo Jo Jo 

The prime on the k integral means that the values kt m+1 l = kM , . . . , kl m l 
must be excluded through a principal-part procedure. We have for the 
interaction of particle j, bearing the charge ej, with the i labeled mode: 



r'3b 
H 



1 



' \ 5 



— e 



.TTTT EE S[ exp{,a[kW. q 0)-2<W]} 
(2?r)2 a=12a=±1 V fcU / 



d ia d 







dpi j) 



.( v (f). e («)M) 2vr 



2^ ] 



(2.15) 



The elements of the metric tensor g have been chosen as g rs = g rs = —5 rs , 
i,r = 1 — * 3. The last matrix element of interest for us describes the 
Coulomb interaction between the two charged particles: 



L 



'(12) 



H 



eie 2 



'dlq(i) _ q (2) 
9q( 2 ) 



-l 











(2.16) 



2.2 Enlargement of dynamics 



We proceed now to an enlargement of dynamics as in previous publications 
fJ, E3|: we multiply the number of variables on physical ground in such a 



way that the original dynamics (2.7) be included as a particular case. The 
choice of a particular enlargement is determined by opportunity linked to 



8 



physical considerations and the properties to be examined |T^]. Since all 
enlargements provide an equivalent alternative description, that degree of 
freedom is welcome. In the present paper, focalized on the self-force on each 
particle, our choice is to define the self-field with respect to each particle. If 
the interest bears on the field far from the two particles, defining the self- 
field with respect to both charged particles would be an alternative useful 
option. 

The elements of the enlarged dynamics will be noted by a supplementary 
upper index tilde "~" , as well for the variables as for the evolution operator. 
Our aim is indeed to eliminate explicit self-interaction processes from the 
evolution, while taking their effect into account. We distinguish formally 
between 5 varieties of oscillators, based on the recognition of self-energy 
parts in the evolution. To each oscillator [i], we associate a discrete index 
that determines which interactions are possible for the oscillator (the index 
j takes the two values 1 and 2). 

[i(sj)] will be the label of an oscillator which has previously interacted with 
the particle j and will further interact with it in a future, without interaction 
with the other particle (j' ^ j), and without playing a role in a measurement: 
by definition, such oscillator does not play a role in the computation of the 
mean values. 

[i(ej)] will be the label of an oscillator which has previously interacted with 
the particle (J) and will no longer interact with it directly: its next inter- 
action should involve the other particle (j 1 ), or it should contribute in the 
computation of mean values. 

[i(f)] will be the label of an oscillator mode which has not previously in- 
teracted with the particles (1) or (2). Its excitation has its origin outside 
the two charges and such an oscillator is free of constraints on its interac- 
tions: either with one of the particle or with an external devise. It provide 
a contribution in the computation of mean values. 

The free evolution of those oscillators is the same as in the original 
dynamics and does not involve a change in their nature. 

The vertices for the computation of < ll[m]|£^|ll[m] > involve Lrj for 
all i : 1 — > m. the numerical value will be preserved for the non-vanishing 
elements. We have to take into account the (possible) change of nature 
of the oscillator after the interaction. We introduce indices corresponding 
to the transition of nature of the field {i{e\f) means that a free oscilla- 
tor i(f) becomes of the emitted e% variety) and we have the non-vanishing 
possibilities: L^ f) \ L^ f) \ L^^\ L^^\ 
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£'l[i(«ea)] l'2[i(s 2 ei)} £'2[i(« a «2)] £'l[i(ei«] |/2[i(e 2S2 )] whUe the 

elements 'l 1[i(eiei] i' 2 ^ 6262 )] j^M^i] £'2[i(- a e a )] '^'i[i(«i« 2 )] £'2[i(*,«i)] 
Lft l ^ eiS2 ^\ £^t*( e2S1 )l vanish by construction. 

Other elements, such as L^ (//1 , L^™, L^^, L^ 1 , 

f'2[i(/e„)] ?'l[i(/«2)] f'2[i(/«i) f'We 2 ] f'2[i(/ei)] f'l[i(e 2 e„] f'2[i(e iei )] ;'l[i(s 2 e 2 ] 
H ' H ' H ' H ' H ' H ' H ' H ' 

;'2[i(siei)] ?'l[i(e 2 ei] ;'2[j(eie 2 )] p'l[i(s 2 ai)] ?'2[i(ais 2 )] ?'l[i(e 2 si)] ?'2[i(eis 2 )] 
H ' H ' H ' H ' H ' H ' H ' 

?'l[i(«2ei] f'2[i(aie a )] f'l[*(s 2 /)] f'2[i(«i/)] f'l[i(e a /)] ?'2[i(ei/)] f'l[i(* 2 s 2 )] 

Z2 [i( * lSl)I , L^ (e2S2)1 , L^ (eiSl)1 vanish obviously since the final label of the 
oscillator does not bear the name of the interacting particle. 

The vertices for the computation of < llfrnJI/^lllfm'] > (m' / m) 
involve a (m + oscillator mode and its disparition from the explicit de- 
scription. The value of the vertices involved, corresponding to L^- 1 is the 

same as the value of L^ m+1 ': we have to consider the non- vanishing pos- 
sibilities for the nature of the (m + oscillator. The oscillator on which 
the integration is performed is considered belonging to the self variety and 
< ll[m]|£#|ll[m'] > will thus involve the following elements: £^ m+1 ( si ^)l j 

£'2[m+l(s 2 /)] £'l[m+l(sie 2 )] £'2[m+l(s 2 ei)] £'l[m+l(sisi)] ^2[m+l(s2S 2 )\ w ^ 

the elements of < ll[m]|£' H |ll[m'] > involving the L^ [m+1(siei)1 and L^ [m+1(s262)1 
vanish by construction. The other oscillators (1 — > m) are unchanged by 
the transition vertex. 

2.3 Constitutive relations-Equivalence conditions 

Matrix elements of the evolution operator for an enlarged dynamics involve 
now the five varieties of oscillators. We have to connect the elements of 
the extended dynamics to the original one. The simplest case involves one 
oscillator only, from the first equation of the hierarchy: 

oo 

dtfuim] = E < n[m]\£ H \n[m'} > f 11[m , } . (2.17) 

m'=0 

For m = 0, we take obviously /n[o] = /n[o]- That first equation means: 

dth m =< 11[0]|£ H |11[0] > / 11[0] + < 11[0]|£ H |11[1(/)] > 
+ < 11[0]|£h|H[1(si)] > Ai[i( Sl )]+ < H[0]|£//|ll[l(s 2 )] > /ii[i( S2 )] 
+ < ll[0]|£tf|ll[l(ei)] > / 11[1(ei)] + < ll[0]|£ ff |H[l(e2)] > /ii[i(e 2 )]- 

(2.18) 
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Ch is composed of the parts part C l H and C 2 H according to the particles 
involved in the interaction. C H acts on /nhmi, /n[i( si )], /ii[i(e 2 )l while C 2 H 
acts on /nfi(/)], /n[i(s 2 )], /ii[i(ei)l- Since we have to recover the equation 

d t f im =< n[0]\C H \ll[m'} > f n[0] + < ll[0]|£ H |ll[m'] > /n[i], (2.19) 

we are led to the constitutive relation Q , (II 

Al[l] = /ll[l(/)] + /ll[l(ei)] + /ll[l(ea)]- ( 2 ' 2 °) 

Indeed, if the conditions /n[i( si )] = /ll[i(ei)] an d /ll[l(s 2 )] = /ll[l(e 3 )] are 
satisfied at the initial time (equivalence conditions), they will remain satis- 
fied for all times and we recover ( |2.19| ) as a particular solution of our set of 
equations. 

Let us consider now the next equations of the hierarchy. 



flt/ll [1(f)] 




< H[l(/) 


\C H \U[l(f)] >/!![!(/)] 






+ 


< H[l(/) 


|£ H |H[2(//)] > /ii [2 (//)] 






+ 


< H[l(/) 


|£ Jf |ll[2(/ Sl )]>/ u[2(/si)] 






+ 


< H[l(/) 


|£ Jf |ll[2(/ S2 )]>/ u[2(/s2)] 






+ 


< H[l(/) 


\C H \ll[2(f ei )} > /ii [2 (/ ei)] 






+ 


< H[l(/) 


|£ H |ll[2(/e 2 )] > / 11[2(/e2)] , 


(2.21) 


<%/ll[l(si)] 




< H[l(si) 


|£ H |ll[l(si)] >/ii[i sl] 






+ 


< H[l(ai) 


|£tf|ll[l(/)] > 






+ 


< H[l(si) 


|£if|ll[l(e 2 )] > /ii[i( e2 )] 






+ 


< ll[l(si) 


|^|ll[2( S i/)]>Ai [a(si/)] 






+ 


< H[l(fli) 


|£//|ll[2(si/)] > /n[ 2 ( SlSl )] 






+ 


< H[l(ai) 


|£h|11[2(s 1 s 2 )] > /ii[ 2 ( SlS2 )] 






+ 


< H[l(si) 


|£ jH -|ll[2(siei)] > /u[ 2 ( siei )] 






+ 


< ll[l(«i) 


|£j3-|ll[2(sie 2 )] > /u[ 2 ( sie2 )], 


(2.22) 


and a similar expression for c?t/ii[i( S2 )] ■ From the equality of the matrix 



elements, we have also (9t/n[i( ei )] = ^/n[i( Sl )] and 9t/ii[i(e 8 )] = 9t/ii[l(« 3 )]- 
Those relations have to be compatible with: 



d t f u[1] =< 11[1]|£ H |11[1] > /ii[i]+ < 11[1]|£h|H[2] > /u [2] . (2.23) 
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We have, for the terms diagonal in the numbers of oscillators: 

<%(/ll[l(/)] +/ll[l(ei)] + fll[l(e 2 )])diag 

=< 11[1(/)]|£ H |11[1(/)] > /ii[i(/)]+ < 11[1(si)]|^h|11[1(si)] > /n[isi] 
+ < 11[1( Si )]|£h|H[1(/)] > /ii[i(/)]+ < ll[l(*i)]|£ff|ll[l(e 2 )] > f 11[1(ea)] 
+ < ims 2 )}\C H \n[l(s 2 )} > / 11[ls2] + < 11[1( S2 )]|£ H |H[1(/)] > /ii[i(/)] 
+ < ll[l( S2 )]|£ H |H[l(ei)] > fume,)] (2-24) 

For these terms that do not involve a field oscillator, that equation is man- 
ifestly compatible with the previous one. Let us consider the other contri- 
butions involving C\j. We have its action on /u[i s i], /unmi, /n[i( e2 )] ' ana - 
this is compatible with the original equation, thanks to the constitutive re- 
lations and to the numerical identification of with /n[i e i] inside the 
equivalence relations. The other terms can be treated in a similar way. 
That relation ( |2.20| ) can be easily generalized for two or more oscillators: 

/ll[2] = /ll[2(//)] + /ll[2(ei/)] + /ll[2(e 2 /)] + Al[2(/ei)] + /ll[2( ex e x )] 

+ /ll[2(e 2 ei)] + /ll[2(/e 2 )] + /ll[2(eie 2 )] + /ll[2(e 2 e 2 )] (2.25) 

and similar expressions for the set of all elements {/iim}- 

The consideration of initial conditions that do not satisfy the equivalence 
conditions requires that the generators of the Poincare Lie algebra be con- 
structed for the extended dynamics and that point is beyond our aim in this 
paper. However, we believe that the distinction between the self-field and 
the external field resists a Lorentz transformation and therefore no problem 
should arise from the extension of dynamics. Moreover, the subdynamics 
operator IT has been proved (in another realisation but within a similar 
framework) by R. Balescu and L. Brenig [|l4| to be relativistically invariant. 
Nevertheless, as far as the adequate construction of the Lie brackets for the 
ten generators of the Lorentz group has not been performed, we have to re- 
strict ourself to the equivalence case. When the compatible (or equivalence) 
conditions are fulfilled, the new dynamics is simply a reformulation of the 
original one. 

2.4 Factorisation properties 

The free fields variables are by definition not connected with the charged par- 
ticles variables. Therefore, the initial conditions concerning the free fields 
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and particles variables can be chosen as independent. The vacuum com- 
ponents {/n[i]} can be factorised into a part, describing the particles and 
the oscillators of the emitted e\ and e2 varieties, and a part describing the 
free variety of the field, for instance an incident field that may be or not 
vanishingly small. We shall therefore write for instance: 

/ll[l(/)] = /ll[0]/[l(/)] 
/ll[2(//)] = /ll[0]/[2(//)] 

/n[2( ei /)] = /ii[i( ei )]/[l(/)] (2-26) 

The free field /[i(/)] distribution function may be in particular the distri- 
bution function fX/f\\ describing the absence of field (field vacuum) and 
considered later on for the computation of the effect of the self-field on the 
motion of the charged particles. In the extended dynamics, the natural 
choice is to consider for fX(f\\ a distribution function corresponding to a 
field of null amplitude and no phase dependance. In those circonstances, 
the function /n[i( Sj )] receive e.g. contributions form {/n[o]/[n(//.../)]} di- 
rectly through the creation operator < ll[l(sj)]|C'|ll[n(// .../)] > (Other 
contributions are written in §4). 



Even in absence of field, the factorisation ( |2,26| ) is not equivalent to a 
factorisation /n[i] = /n[0]/|i(/)]i n the original dynamics. In the equivalence 
conditions, we have numerically that the functions /n[i( Sj )] and /n[i( ej )] 
coincide. The constitutive relation ( 2.20|) requires that /n[i] = /n[i(/)i + 



/ll[l(ei)l + /ll[l(e 2 )l- Therefore, in the original representation, we are not 
allowed to consider a factorisation /n[i] = /n[o]/[Y(/)i corresponding to the 
absence of field at the time considered. If we impose at some time /n[i] = 

/ll[0] /[!(/)]> at tne same time > we nave to consider /n[l(/)] = /n[o]/[Y(/)] ~ 
/n[i(ei)] ~~ /n[i(e 2 )]- Therefore, we have to admit the presence in fimrfsi of 
contributions — /nfl(ei)] — /n[i(e 2 )]- When computing the equation of evolu- 
tion of the charged particles, those terms play a role directly through, for in- 
stance, the element < 11[0]|£# |11[1(/)] > of the first equation of the hierar- 
chy ( |2.18| ) . That contribution has to be combined with the contribution aris- 
ing from the kinetic operator ©. Upon that imposition /n[i] = /il[0]/[Y(/)]> 
it is mandatory to consider also that contribution to have a valid comparison 
with the usual results. 

In conclusion, we have been able to provide a set of evolution equations 
for the reduced distribution functions that enables the explicit identification 
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of the self-field. 



3 The kinetic operator 

3.1 General concept 

We now proceed to the construction of a subynamics associated with the 
enlarged dynamics. First of all, we have to define the vacuum and correlation 
states. A correlation state contains at least one self oscillator while the 
vacuum (of correlation) is defined as the set {/nm} where all oscillators 
are of the free / and emitted e\ and ei varieties. The construction of the 
subdynamics rests on that distinction. All the formal results of the Brussels 
group, concerning its construction rules and its formal properties, are applied 
directly, with our specific realisation of the operators involved. 

Inside the subdynamics, the vacuum components obey close evolution 
equations, namely kinetic equations. These equations are not time reversal 
invariant, hence their name. For their determination, we can limit ourselves 
to the consideration of the vacuum-vacuum elements of the superoperator 
£(t) 0, its t = value defining the IT operator. We take for granted the 
usual properties of idempotency, factorised structure and commutation of II 
with the evolution operator C. The explicit verification of those properties 
requires the explicit knowledge of all the elements of £(£) Jl3| , and not only 
of the vacuum- vacuum ones. 

In the enlarged dynamics, the evolution equation takes the form: 

oo 

dj n[m ] = E < 11H|£h|HK] > / 11K] , (3.1) 

m'=0 

which involves all the varieties of the oscillators. The kinetic operator 
associated with the subdynamics provides in an exact way close equations 
involving the vacuum oscillators only: 

oo 

Wiim = E < iiMI®|uK] > v&iw ( 3 - 2 ) 

m'=0 

The value of that operator can be reached by the direct computation of the 
vacuum-vacuum elements of the superoperator X(t). 

The hierarchical form of the equations (3.1) (m' > m) enables the deter- 



mination of the elements of in a successive way. The elements of G that 
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do not involve an oscillator are the same as those of jCh, and therefore also 
the same as 

We proceed to the computation of the first non diagonal element 

< 11[0]|G|11[1(/)] >. It is based on the evaluation of the corresponding ele- 
ment < 11[0]|S(£)|11[1(/)] >. Its evaluation is performed in a perturbative 
way. The elements will be affected by a couple of upper indices which de- 
scribes the number of Coulomb interaction and the power of interaction with 
the oscillators. The simplest element is of course < ll[0]|S(i)|ll[l(/)] >( 0,1 \ 
in which the Coulomb interaction is not considered and only one interaction 
with the (free) oscillator takes place. Such element involves no self oscillator 
and we have trivially: 

< ll[0]|E(t)|ll[l(/)] >(°' 1 )=< ll[0]|exp£ H t|ll[l(/)] >(°' 1 ) . (3.3) 

The vacuum- vacuum elements of 2(0) are noted A and from the general 
relation valid for vacuum- vacuum elements 

< ll[0]|E(t)|ll[l(/)] >=< ll[0]|e § U|ll[l(/)] >, (3.4) 

we have < 11[0]|6|11[1(/)] >W=< ll[0]|i H |ll[l(/)] >. 

3.2 Construction rules of Vt{t)V 

The general construction rules are illustrated on the first non trivial element 

< 11[0]|E|11[1(/)] >( ' 2 ), in which the Coulomb interaction is not considered 
and two interactions with an oscillator take place. Such element involves 
one self oscillator if the two interactions involve the same particle. If they 
involve different particles, only physical states are present in the contribution 
and we have anew the equivalence of the corresponding elements of and 
Ch- We dispense ourself of a supplementary index and concentrate on the 
contribution involving a self oscillator. We have: 

< 11[0]|E|11[1(/)] >(°< 2 )= -± fdze-** £ 

AmJ ° 3=1,2 \*-*£jff/ii[0],ll[0] 

x i < ll[0]|£ff|ll[l(sj)] > 



\ z ii[i( Sj -)],n[i(«i)] 
x z<ll[l( Sj )]|^|ll[l(/)]> ■ (3-5) 
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The prime on the integral sign means that only poles corresponding to prop- 
agators arising to vacuum states (without self oscillators) have to be included 
in the path c. In the present case, the pole due to the intermediate prop- 
agator is thus excluded from the path. That selection of poles corresponds 
to the recipe to construct the subdynamics. The accidental coincidence of 
poles due to the correlation and vacuum propagators is avoided by adding 
a positive imaginary infinitesimal ie to the correlation propagators when 
computing the residues. Another formulation of the recipe is the following: 
a positive imaginary infinitesimal ie is first added to all propagators corre- 
sponding to the correlation states and the path c encloses then the real axis, 
above — ie. When no resonance can occur, the ie can be dropped. 

3.3 Fourier representation 

The evaluation is more easy in variables such that the free motion opera- 
tor is diagonal. For the free motion of particles, those variables are well 
known and correspond to the Fourier transform of the original spatial vari- 
ables. Therefore, we will replace the unknown /n[o] where the variables 
xW, are (q^,p^), (q*- 2 -*, P^) by new functions depending on variables 
(k^pW),^ 2 )^ 2 )). We will not introduce a new symbol: the nature 
of the argument specifies the function under consideration. The transition 
between the two functions is provided by (we use Balescu's choice for the 
normalisation factor): 

/ii[0](q (1) ,P (1) ,q {2) ,P (2) ) 

= JL/ ^^^^e^-^^^^n^^kW,^),^)^^)), 
/npiCkW.pW.kM.pP)) 

= ld 3 qi d 3 q2 e^ (1) ^ (2) ^ (3.6) 

All functions fu\ m ] have to be similarly replaced. 

We have to perform a similar change with respect to the variables as- 
sociated with the oscillators. As the functions are periodic in the variables 
£M ; Fourier series are relevant. The function /n[i] becomes a new function 

depending for the oscillator on the new variables (r]a\rria\a = 1,2) (m^ 
discrete) in place of the continuous variables (rja\ £a\ = 1,2). 

/nmCkW.pW.kM.pWs^UfU™,^ 1 ;^ 1 ) 
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x /n W (k« k( 2 ) , p( 2 ) ; ^U! 11 , 4 1] , tf 1 ; kW ) , 

= e-^rffN"^ 11 ) 
[i] [i] 

mi L J it) L J 



(3.7) 



x / n[1] (kW , pW , k( 2 ) , P ( 2 ) ; r,? , m? , 4 ] , ; k N ) , (3.8) 

where the summations run on all integers, positive and negative. 

In Fourier variables, the one particle and one oscillator free motion op- 
erators take a simple diagonal form: 



L 



o(i) 

H 



2 



while we have for L 



'(12) 
H 



H 8vr 2 7 

or (alternative more usual form) 



a 



d 



dpW dpi 2 ) 



d 



d 



d P w dpi 2 ) 



We have for L% [i] : 



(3.9) 



(3.10) 



(3.11) 



L 



'M 

H 



1 



_ \ " \ " / 



dpi 



9 a 



M-«) 



x exp a 



-kW. 



9k0) dm l 



[i] 



(3.12) 
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The only difference is the replacement of the variable qf by the par- 
tial derivative — i— 7-y and a similar transposition for the angle variable of 

dkf. ' 

the field. The notation exp —a \ i enables to take into account the non- 

drn a 

diagonality of Ljj 1 with respect to the index m a : the transition is ±1 
according to the value of a. Another possibility is the introduction of the 
factor V i H] S m , writing with a prime the corresponding argument 

m l a ' ,nia '+a 

of the function on which the matrix element acts. 



3.4 First order kinetic operator 

In those variables, the operator < 1 1 [0] 1 1 1 1 [1 (/)] >(°' 1 )=< 11[0]|£h |11[1(/)] > 

takes a simple form, due to the presence of a front factor 5 m „<5 [11 „, a 2 = 1 

m\,0 m 2 ,0 



2— — 



and the property 

•/ 

when acting on a regular function. 

< 11[0]|G|11[1(/)] >(°' 1 )= - 2 ej r / d 3 fcW / drjp / cfy 

jt^ 2 (2vr)2 •/ ./o 

E ^l/m^O E E p 
[i] m Ii] 1 2 «=1,2«=±1 \ K / 



.[1] 



X 



5 



dp 



(j) 



x exp a ■ 



d 



(3.13) 



3.5 Recovering the Lorentz force 



For pedagogical reasons, we show in Appendix A that under conditions 
of independence (factorisation) of the field and one particle distribution 
function, describing a particle sharply located at r(t), this expression leads 
to 



dj(k,p,t) i = -e < E ± (r(t)) +p x B-^t)) > .V p /(k,p,t). (3.14) 
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e < E 1 (r(t))+vxB 1 (r(t)) > is the usual electromagnetic force acting on the 
particle. The minus sign is easily accounted for. If the distribution function 
corresponds to a well defined value of the velocity and an uniform acceler- 
ation a, it can be written /(k,v, t) oc <5(v — vo — at) and for the time de- 
pendence due to the acceleration, we have dtf(k, v,i) = — a.V v /(k, v, t). 
That expression clearly shows that the present formalism yields the well 
known expression for the Lorentz force. That relation ( |3.14 ) will be used in 
§ 6. 

3.6 Evaluation of the second order < 11[0]|E|11[1(/)] > 

While the evolution of the field will be considered in the next section, we 
focus here on the reaction of the particles to the presence of the field, namely 
the radiative corrections to the direct interaction between the particles and 
the field. To provide their contribution to the evolution of the two-particle 
distribution function, only the elements of the kinetic operator that acts 
on the (reduced) distribution function of the two charged particles and one 
mode of the field are relevant. When acting on the distribution function 
corresponding to the absence of field (defined in the extended dynamics), 
they will determine the first radiative correction to the free motion of the 
particles. 

The operator < 11[0]|E|11[1(/)] > (0 ' 2) Q in the new variables is now 
determined. The two interactions have to involve the same particle. We use 
the same convention for the index of particles as for the polarisation of the 
oscillators: f is 2 when j is 1 and vice versa. We get: 

< 11[0]|E|11[1(/)] > (0 ' 2) 

= 2riL dZ€ ~ lZt ^ 2 C-kW.v(;)-k(i'). v 0'')) ( "' )6j ^| 

/p OO P oo 

«"| *F E E 

_M_[l] a=l,2o=± 



x i. 



x cxp a 



=±i 

d 

dpi j) 

d d 1 



kW. 



19 



kO).vO') - kO v ). v O') + m(m l a ] + 



1 \ 2 



x(— i)ej 



0=1,2 6=±1 



E E ( fc fl] 



-7T 



(v y) .e^ )[1] ) 



dp® 



(roj? - b) 



exp&i -kW. 



d 



d 



kO').v(i) - kOO.vO'') + fcWrf +mL 1 , 1 ; 



(3.15) 



The displacement operators can be transferee! at the right of the expression 
to provide 



< 11[0]|S|11[1(/)] > 



(0,2) 



(2vr)t 



= — f dze~ lzt V (z-k^.v^-k^.v^'A \-i)e j — 

//* oo /* oo 

4 1] / E Vio^io E E 

a=l,2o=± 



km 



=±i 







( 2 _ k o). v o) + okW.v^') - k^").v^") + fciVL 11 + 4 1 ' - «)) 1 



1 



[1]\ 5 



x(-i) ej — E E (rm 



- - "V (e t , K r , -e r , « t , ;j 



_^ (v 0'). e (/3)[l]) 



/ _d_ b_ 



1 

V dr >/3 Vf3 



(z - k^.v^ + &k(M).vk') + okW.v 
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_ k (i') . v 0") + k W ( m W + m Pj _ 6 _ a )) 



exp < -(a + 6)kW. 







a 



drn^a 



(3.16) 



The summation over m 



i i 



and the Kronecker delta functions of 
the variables = and mjjj = are also written at the right of the 



i' 1 



expression. In the explicit computation, a separation has to be performed 
between the contributions with (3 = a and (3 ^ a on one hand, a = b and 
a 6 on the other hand. In our future expressions, the first sign = or ^ will 
refer to the polarisation index while the second one to the relative value of 
a and b. 

From that expression, the computation of < 11[0]|S|11[1(/)] >(°> 2 ) is 
rather straightforward: we have to proceed formally to the derivatives with 
respect to the mechanical momentums of the particle. A further index can 
be introduced to reflect the factor on which they act. The last operation 
is an integral by residue that takes into account the poles due to the first 
and last propagators that correspond to the vacuum of correlation. For the 
sake of illustration, let us consider one of the contributions that will matter 
when acting in absence of free field: 



< n[0]|E|ii[i(/)] >i°; 3 % 







E E 

a=l,2a=±l 




1 



kCtf.vC?) + akW.vC?) - kCiO.vC*"') - afcW 




1 



k0'). v 0) -kO'). v (j') 







dp 



Cj) 



»[1] ,,(?) />)[!] L.W 9 
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z-kO').vO') -kO'O.vO'') 



[jfe[i] e WW-/*' t ,0-) ( 4«)[i] Jb W_ e W[i]^]; 



9 



<9p 



(J) / 



or 



[!] I 1 ] 



(3.17) 



In that expression, the partial derivative — yjy acts on everything at its 



right: the factors v~! and (in the first term), the factors v s , and 



dv 



— fiy I in the second term and, for both terms, the momentum depen- 

9p y r v ) 

dence of the distribution function on which < 11[0]|E|11[1(/)] >=3^ is 
applied. The only singularity to be included in the close path c is the 

7y ) . No coincidence is possible with the pole of 



pole of ^ 



i 



2 -k(3').v(j)-k0').v0'') 
1 



• i i n\ — ftn — tmt (v^ < 1) and no ie has to be introduced 

here. The obtention of < 11[0]|£|11[1(/)] >2?3^£ is therefore straightfor- 
ward. 



< 11[0]|S|11[1(/)] > {0,2) 



3,^ 

i 



a=l,2o=±l V K 



(-it) 



[1] 
: 2 



akW.vO) - oifeW 



1 



-akW.vO') +afcW 



5 



9 



- .(v^leWW) 2 



+ 



1 



akW.vC?) - akW 
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2(-it) 



akW.vOO -ofcW 



+ 2 



1 



9 



a- ; 



[i] [i] 



(3.18) 



3.7 The second order evolution operator 



From the general properties of the sub dynamics, we have the links (3.4) 
between the subdynamics operator S(t) and the evolution operator G for 
the vacuum-vacuum elements. Since we have a limited choice of possibilities 
for the vacuum intermediate states, we have 



^<11[0]|E|11[1(/)] >L°3 2 ^ 



t=0 



= < 11[0]|6|11[1(/)] >L°3 2) ^< H[l(/)]|i|H[l(/)] > ( °' 0) 

+ < 11[0]|£°|11[0] X 11[0]|.A|11[1(/)] • (3.19) 



Using < 11[1(/)]|^|11[1(/)] >(°.°)= land< 11[0]|£°|11[0] >= -i\k.W.v®+ 

),2) 



k0"). v 0")] 5 we can obtain directly < 11[0] |0| 11 [!(/)] >^' 2) 



< ii[o]|e|ii[i(/)] > {0,2) 



3,^ 



x (-0 



1 

akW.vO') -aifcM 
d 



[fcWe? [1] -/%5 ) (i a)[1] fc; i -e?%f 1 )] 



9 



dp9 



-7r(v^.eWW) 2 



1 



okW.vO') -afcW 
9 



Op 



U) 
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•SI 



(3.20) 



[i] [i] 



All contributions are treated in a similar way and recombined such that the 
conservation of the norm is manifest (due to the front factor ^ ): 



< 11[O]|0|11[1(/)] >tf 



= E 4t [<p*p r dv ? f4 1] E iff 



X 



akW.vCi) -akin). 



a=l,2a=±l 

fcW4 a)[11 -A< 3) (ef )[1, 4 1 ]- e w[ 1 ifci 11 



E ^1,0^,0 



[1] [1] 



1 5 



1 ^ 2 



akW.vO') - oifeW 



E ^w.o^w.o- 

rn rn 1 z 



(3.21) 



,m 2 



Here, the commutation property of [feMe^*^ 1 ' — g st v^\e[ a ^ k^ — e^^fe} 1 ')] 



and — ttt has been used 

dpi 11 



OPr 

For the sake of completeness, similar expressions are provided in ap- 
pendix B for < 11[O]|0|11[1(/)] >L°4 2) , < 11[O]|0|11[1(/)] >^L 2) and 

< ii[o]|6|ii[i(/)] - 



(0,2) 
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3.8 Evolution in the absence of incident field 



The knowledge of the operator < 11[O]|0|11[1(/)] >(°' 2 ) enables us to look 
for the behaviour of the particles when evolving into the absence of exter- 
nal field. Other matrix elements could also be considered to provide, for 
instance, the vertex (charge) renormalisation due to the self-field but are 
outside our scope in this paper. 



We can use the factorisation ( 2.18| ): 

Al[l(/)] = Al[i 



]/[!(/)]' 



(3.22) 



and use for /[i(/)] the distribution function /rw/)] corresponding to the ab- 

and t]2 



sence of field: namely the limit for rji 



rji)6(r). 



[i] 



-> of: 



Kr 
i] 
i 



m. 



mk J ,0 



(3.23) 



Since the variables r/[ and 7/2 are integrated from to 00, the limit r\\ — > 
and r\2 — > has to be taken after that we have performed that integration. 
We identify in < 11[O]|0|11[1(/)] >(°' 2 ) all the terms which could provide a 
non vanishing contribution. The summations over mj 1 and provide a 
vanishing result if a displacement operator on and is involved: we 
would then meet a product of Kronecker's delta functions with incompatible 
arguments. Therefore, the only possible non-vanishing ones would arise from 
the contribution < 11[0]|G|11[1(/)] Since that contribution involves 

a front factor r/a\ the presence of 5(rj^ — r]i)S(r]^ — %) for 771 — * and 
??2 — ► provides a vanishing result except for the contributions in which the 
derivative of the Dirac's delta function appears. The second term in ( 3.2l|) ) 
provides therefore a vanishing result and we are left with 

(0,2) 



< 11 



=1,2 



l©lll[l(/)] 

3 (2vr) 3 dp (?) 
1 



>_ 



■f [!(/)] 



d 3 kW / 
JO 





Kr jrKr 

r- 
x ,0 mi. 



(3.24) 
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The summations and integrations over the fields variables can be performed 
in a staightforward way using the Kronecker's and Dirac's delta functions 
(after an integration by parts and performing the limits 771 — > 0, 772 — ► 0) 
and we have: 

< ll[0]|O|ll[l(/)] >L% 2) /J (/)] 

= ^ 2 H)e W^F y / d3fc[11 ^2aSi W («kW.vO')-afcm) 
x [ fc [i] e (-)[i] _^^) (e WW4i] -eWWfcfl)] [_27r(vW. e< a >W)] . (3.25) 

Nevertheless, this last expression vanishes also by parity for the summation 
over a. 

3.9 A useful expression :< 11[0]|.A|11[1(/)] >(°' 2 ) 

For future use, the non vanishing element of < 11[0]|^|11[1(/)] > (0 ' 2) /^ (/)] 
is required: 

< 11[0]L4|11[1(/)] >(°' 2 ) fj (f)] 



3 

X 



<9 



The derivative with respect to pi- acts of course on all possible dependences 
at its right. 

3.10 Physical interpretation 

The concept of renormalised mass has now to be extracted from the kinetic 
equation, by combining radiative corrections with the free motion operator. 
Our expression of the second order element < ll[0]|O|ll[l(/)] >tf of the 
kinetic operator shows that it vanishes for a free particle that is not accel- 
erated by external fields nor a Coulomb interaction No mass correction is 
provided. Indeed, if a relativistic expression is used for the energy of the 
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particle, the momentum and energy conservations cannot be simultaneously 
satisfied by an emission act of a non-accelerated charged particle. No reso- 
nant process is possible and their absence implies that no ie is required and 
the propagor is odd in a. 

This result could be expected from general considerations from our knowl- 
edge of general properties of the subdynamics || when the propagator in- 
volved cannot be resonant. In Brussels terminology H, the second order 
kinetic operator vanishes for parity reasons: it is well known that, at that 
second order, the contribution to the kinetic operator (called ^2) arises from 
a Dirac delta "function" and not from the principal part of the (usually reg- 
ularised by a ie) propagator. No regularisation has been required here when 
acting in absence of field and the kinetic operator provides a vanishing con- 
tribution. 

The effect of the coupling with the field vacuum is therefore to be 
searched in other terms. Indeed, radiation emission is present when the par- 
ticles are accelerated. We consider as afirst step in §5 < 11[0]|6|11[1(/)] >( 1,2 ) 
The acceleration provided by the Coulomb interaction will induce a back re- 
action on the motion of the particles. 

The next section §4 is devoted to the computation of the field generated 
by the particle in their free motion. 



4 Emitted field in a free motion 



In order to get a better understanding of our previous result ( gjgp , we intend 



to analyse the field emitted by the particles and their back reaction at the 
same order (0,2) in the interactions. The contribution is the same as for two 
particles moving independently. We will use the equivalence conditions that 
enable to get the emitted field from the self-field. The subdynamics theory 
determines the last one through the so-called creation operator C. Since we 
know that S(i) = Ce &t A for the correlation- vacuum elements, we focus on 
the elements of £(i)( 0,1 ) that provides a contribution to 

/li[i(isj)] =< ll[l(*i)]|Clll[l(/)] > 

+ <H[l( ai )]|(7|ll[l(e/)]>/ 11[W)] 

+ < ll[l( ai )]|(7|ll[2(//)] > f n[2iff)] + . • • (4.1) 

Since numerically, in the equivalence conditions, we have the equality fn[l(ie )] = 
/n[i(i Sj .)], the lowest order contribution to /n[i( iSj .)] requires < ll[l(sj)]|£(i)|ll[l(/)] > (0,1) 
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that determines the lowest order contribution to the creation operator. Only 
the terms that provide a contribution when acting in absence of field are 
considered. 



4.1 Correlation- vacuum element of 

In place of ( p,15| ), we start from 



< ll[l(^-)]|S(t)|ll[l(/)] ><o,U jv 



i(/)] 



2m 



x — i e 



dz e 



1 



-izt 



z- kO').v(i) - kO'O.vO'O + fcW(mL J + m^ J ) 



J (2vr)l,§ 2a S 1 l^ 



[1]\ 2 



dp® 



-7rfv^).e^W 



: o [1] [l]( m /3 ] °) 
1 



exp a 



-kW. 



9 



9 



z - kO').v(i) - kO'O.vCf') + fc'^mj 1 + mJJ) 



J 1 ! 



[l] 





(4.2) 



where we have taken into account that the final state ( < ll[l(sj)]|) is 

a field correlated state, hence the presence of feW(ma + ^45) i n the first 
propagator. Only the pole of the second propagator (due to a vacuum state) 
is enclosed by the path c. Moving the displacement operators to the right 
and using afterwards + 



my + mpl = m a + mj = a, we get 



< ll[l( Si )]|E(t)|ll[l(/)] >(°' 1 ) fj (f)] 
-1 



2vri 
1 

(2tt)! 



dz e 



-izt 



E E 

/?=l,2a=±l 



1 



kCjJ.vCjJ-kCjO.vO'O+oJfeW 



2 • 



-i 



3pC>) 
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-Tr^.e^W) 2 



d 



kZ _ ( k (j) _ O k[i]). v 0") - kOO.vO'); 

:eXp " ( akI11 w) ^ - ^? ] - ^)^U^,o- ^ 43 ) 

This expression can be computed easily and is identified with the same order 
of<ll[l( Sj )]|Ce0*i|ll[l(/)]>/f (/)] . 

4.2 First order creation operator 

Therefore, 

< ll[l( Sj )]|(7|ll[l(/)] >(°' 1 ) fj {f)] 

1 



= E E 

/9=l,2o=±l 



-akW.vC?) +akW 



/ I 1 ! \ 5 
1 I V f3 



-7r(v^.eWW) 2 



5 



exp — ak 



"ak(?) 



[i] 



+ E E 



-i 



/3= t 2a fe(-« k[11 - v(j) +^ [11 ) 2( ' )ej (2vr)l ^ 



1 / ^ M 



9;; 



(i) 

's 



dp 1 ', 



x «"> - ( okl11 w) J <"™ " *> i <"™ " <«> 

The limits rip — > 0, 77^ — > have to be performed after the integration over 
r$ and r^ 1 ,' . 

Using (3.6), the original variables £W are reintroduced in place of mW. 
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4.3 Field associated with free particles 

Since /n[i( e .)] = /ii[i(s 3 -)]> we have at first order in the field interaction and 
zero*' 1 order in the Coulomb interaction: 



/ 



(0,1) 
n[i(e 3 )] 



E E 



-27ri«™) 



(3=1,2 a=±l 



1 



-okW.vO') + aJfcW 



-z e 



J (2vr)l 



9?^ 



-7T 



( V (i). e (/3)W) 2 



9 



exp — a 



x^J} 1 -iip)8(ij$ - ?7/3')/ii[o] 



+ E E* 

/3=l,2a=±l 



-27rt(a{W)_ 



-T 5 



(-akW.vO') +a/fcW) 2 ' 3 (27r)i U [1] 



<9t> 



(./) 



x exp — a 







0kO') 



[i] 



(4.5) 



Combining the form (2.2) for the observables associated with the transverse 
electric field and the form (2.5), we get the average transverse component 
of the electric field and we can write 



< E^(x) > 



e 3 (0,l) 



Jo Jo Jo Jo J 

x-^r E E fc wi e?(kW)^(kW) 
(27TJ2 a =l,2a'=±l 

x exp{ia'[kW .x - 2< Q (k™ )}}f^ ej] (* (1) , ; X [1] ; k [1] ). (4.6) 

We now take for /n[o] in ( fO] ) a distribution function corresponding to sharp 
values of the positions and momenta and (3.6) enables to get the expression 
into the variables k.^\ p^: 

/ll[0](q (1) ,P (1) ,q (2) ,P (2) ) = 5(q {1) -qi)5(q {2) -q 2 )5( P (1) -Pl)5(p (2) -P2), 



30 



/ u[0] (k« , p« , k ( 2 ) , P (2) ) = / A (2) e~^ W -<i (1)+k(2) ^ (2) ) 

x<5(q« - qi)5(q( 2 ) - q 2 )<5(p« - P i)<5(p (2) - p 2 ) 

= e -(k (1) -qi+k( 2 ).q 2 ) 5(p (l) _ pi)(5(p (2) _ p2 ). (4 . 7) 

Combining the previous terms and performing all trivial integrations, we 
get: 

< E^(x) >«M = H )( 2 „) e .JL_ / dS k W £ ^e-(kW) 

1 



exp{ia[ k W.(x — qj)]} 



ikW.v, 



a=l,2a=±l 
,(a)[l]^ 



(4.8) 



The value < E^-(x) > e i(°' 1 ) is now determined by an expression that involves 
the values of the position and momentum pj of the charged particle j at 
the same time. 

The summations over a and over the polarisation vectors lead to 



< E ± (x) > e ^ = H)(27r) ej 
-expl-flkW.tx-q,)]}" ' 



1 



(2tt) 3 
1 



d 3 kW [exp{i[kW.(x-qj)]} 
\ / (vj.kWjkWN 



kW.Vj-fcW 



.(4.9) 



This expression behaves obviously as 



|x-q 3 f 



and does not describe a propa- 



gating field. It presents a discontinuity and vanishes exactly at the location 
of the particle x = q^ since the integrand is identically null for that value. 
In view of the value provided by fl3.14|) , this explains why the corresponding 
terms in the kinetic operator vanishes. In our approach, the contribution 
of each mode to the self- interaction of the electron is computed first. When 
no acceleration mechanism is provided, each contribution vanishes exactly. 
In the usual approach, the emitted fields are computed from the Lienard- 
Wiechert potential and evaluated, via the Abraham-Lorentz model, at the 
localisation of the electron. 

The expression for the electric field can be further analysed. We can write 

Q as the sum of two contributions by replacing ( (yj — - "^[i] 

£H b y 



the sum 



( v i _ fePl) + (ffif 



kli] (v 3 -.kM)kM 



The contribution of the second 
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term is: 



<E ± (x)>^ (0 ' 1) =-(4vr) e ' 



d 3 kW sin[kW.(x- qj )] 



J (2vr)3 

(4 " )ej (2^)3 Vx / ^ [1] cos[kW.(x- q ,)]-i-. 



(4.10) 



From the known relations J d 3 kW exp{— i[kM .x]} = 47rJ? an d 
/ <i 3 kW exp{z[kM.x]} ^ fc ^ 2 = 2-7r 2 i, the following identification is possible: 



<E ± (X) >f j(0 ' 1) : 



ej V x - 



1 



x 



Oil 



< E"(x) > 



(4.11) 



Therefore, the first contribution < E _L (x) >a J 0,1 should be identified with 
the complete electric field < E(x) > e i(°' 1 ): 



ej (0,l). 



(-i)(2*)e r 



1 



<E-(x)> a - v . A ...^ N , 

-exp{-i[kW.(x - qj )]}] (k [1] Vj - kl 1 
1 



d 3 kW [exp{i[kW.(x- qi )]} 
1 



(-i)(27r) ei 



(2*0* 



d 3 kW exp{i[kW.(x- qj )]} 



exp{-i[kW.(x- qi )]} 



1 



fcW -kW.v, 



kW-fcWv, 



(4.12) 



If the x axis is placed along (x — qj) and the y axis along v±j, defined by 
vj_j = Vj — ^'^x^.jf we show in Appendix C that: 



< E x (x) 



2 ^ |x- qj | 2 



(4.13) 



On the other hand (11.154) of [[?J gives us the field in terms of the position 
of the charges at the same time: (r = |x — <jj-|, /3 can be identified with Vj, 
q = ej, cosip = n.^ 2 -, r = rn, 7 2 — 1 • 



E 



qr 



r 3 7 2 (l — P 2 sin 2 ip) 



(4.14) 



Therefore, our expressions ( |4.9|44.12 ) reproduce correctly usual results for 
the complete and transverse electric field outside the location of the charged 
particle. 
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4.4 Physical interpretation 



The field that has just been computed is the field generated by a particle in 
uniform motion since the Coulomb interaction between the charged particles 
or an outside field is not taken into account in its contribution. That field is 
therefore equivalent to the field that can be deduced fom the static Coulomb 
field through a Lorentz transformation and this corresponds indeed to our 
result. This point has to be viewed as a confirmation of the correctness of 
our alternative approach. The usual (relativistic) expression is recovered 
outside the location of the particle x = q^. 

On physical grounds, a charged particle in free motion does not emit a 
field and should experiment no self-force. Our expression ( |4.9| ) is in accor- 
dance with that property since the self-field vanishes exactly at the location 
of the point charged particle. That property holds as well for the transverse 
field < E J -(x) > as for the complete one. This explains why the correspond- 
ing terms in the kinetic operator O ( |3. 25 ) vanish. 



5 The radiative reaction force due to the Coulomb 
interaction 

In order to get a contribution to the reactive force due to the self-interaction 
of the particles, a mechanism of acceleration of the particles has to be pro- 
vided. We have chosen to consider the Coulomb interaction between the 
charged particles as responsible for the acceleration. Other mechanisms are 
possible, such as the presence of a non-vanishing free field, or the consid- 
eration of the field emitted by the other particles but they are not treated 
here. We have to evaluate the elements of E(i)^ 1 ' 2 ^ (corresponding to one 
Coulomb interaction and two interactions with the transverse fields) that 
provide a contribution to < 11[O]|0|11[1(/)] >^' 2 \ when that operator acts 
on the field vacuum. The Coulomb interaction between the two particles 
can occur as the first, the second or the last interaction. Since we know 
that < ll[0]|O|ll[l(/)] >( ' 2 ) provides a vanishing result when acting on 
the field vacuum, we expect that the only non vanishing contribution arises 
when the Coulomb interaction takes place between or after the interaction 
of the particles with the transverse field. 
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5.1 The subdynamics operator 

Therefore, we first focus on (the lower index FPF describes the order of the 
interactions): 

< 11[0]|E|11[1(/)] >g& 

x (—i)ej 



.ff o U-k0-).v0)-k(i'). v 0-'), 

J— 

(2k) 2 J JO Jo m ji, a=l,2a=±l 



5 



x expa < — k' J. 



[i] J? ^_ 



X 



v z - k(i).v(i) - kO'O.vO'O + fc[i](mL 1] + mJJ 



Xejej '8^J dl p l [d^)~d^U 

1 



9 a 



x 



v z - k(i).v(i) - kO'O.vO'') + fc[i](mL 1] + mJJ), 



x (-^;ttt E E (iW 



(2vr)2 (9=1) 2 6=±i 



x 



-vr(v(^e(^])(2-A I -4 T ( f nW-6) 

^ 



9 



exp ^- k[1] W"^ 



(5.1) 

>(0,2) 



This expression is very similar to the expression of < 11[0]|S|11[1(/)] 
( |3.16| ), but with the supplementary factors due to the Coulomb interac- 
tion: the matrix element (3.10) and a propagator. The order of all the 
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elements has to be strictly respected, on view of the presence of displace- 
ment and derivation operators. The contributions due to the different or- 
ders of interaction {FFP and PFF) are evaluated from similar expres- 
sions. We can proceed exactly as for the second order contribution < 
11[0]|E|11[1(/)] >(°' 2 ). A difference is the presence of a denominator that 
can be resonant. The subdynamics theory has prescribed, from the be- 
ginning of its elaboration, that a propagator corresponding to a correla- 
tion state has to be treated with an ie. A second difference is the con- 
sideration of < 11[0]|.A|11[1(/)] >(°- 2 ) fX {f)] (|3T26|) for the extraction of 

< 11[0]|6|11[1(/)] > from < ll[0]|£(i)|ll[l(/)] > f[ 1(f)[ through 

( |3.4| ). Moreover, we have to take into account, in the final simplifications, 
that the matrix elements associated with the Coulomb interaction and the 
field interaction do not commute. The terms corresponding to that case are 
affected by a lower index II, the other ones by an index I. 



5.2 Kinetic operator 



Straightforward but very lenghty computations lead to the following expres- 
sion (as expected, the PFF order of interaction has provided a vanishing 
result) 

< 11[0]|G|11[1(/)] >j' 2) fj m 



E( 

j=l,2 
1 



-i)e^jej>- 



1 1 



I fulfil £ £ 

J J « _ -jo L 



k 



l 



(2vr) 3 4vr 

[fc [l]4a)[l]_ / t u (i) (e W[i]^[ 1 ]_ e ( a ) [ i] fc [l] )] 



a=l,2a=±l 




dp 



U) 



ie + akN.vCi) - ofeW J We + (il + okWj.vO") - il.v(*') - afcW 



xl v (-akW) 



or 



Cj) 



27r(v (i) .e (Q)[1] )e 2 'V 8kU) akU 7 ) J } 



(5.2) 



< 11[0]|9|11[1(/)] >& 2) /7 (/)] 



E 

1 



1 1 

(2vr) 3 4^ 



/ / <* 3 4 E E 



ci=l,2a=±l 



A; 



l 



[ k m e (a)[i] _ g st v ti) {e (<x)[i] k li] -eWWjfej 1 ])]. 



9 



Op 



Cj) 
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l v ie + (il + akW).vO) - il.v(i') - akW J U + akW.vO') - akK J 

fbffl 1 ( 9 9 \ 

xi^ e («)[i] e H^ ikO 7 );. (5.3) 

dpv 

The order of integration respects the ordering of the apparition of the ver- 
tices: in all remaining contributions, the integration over the field modes 
has to be performed after the integration over the wave number 1 exchanged 
by the Coulomb interaction. An opposite order of integration would have 
been required in the contribution involving the order PFF of the vertices. 

For the sake of completion, let us compute (^xjy)- We have p = 

from which we deduce ( ^fl) = 5 -& r 

(l-v 2 )? \dp { v J> J ( m 2+(pO'))2)^ (m 2 .+( p U))2p 

The partial derivative d (i) can be placed in front of the matrix element 

dp)- ' 

(i) 

since the contribution when it acts on the factor vl provide a vanishing 
result. The property can be checked explicitly. Taking into account the 
value of the matrix tensor g st = —5^[, the first term of the derivative, 
with a , can be seen to involve the scalar product of the vectors e( Q )M 
and kM that vanishes by definition of the polarisation vector. The second 
contribution, with the product Pu Pv , vanishes by symmetry. This result 
is not unexpected and reflects that the magnetic force is orthogonal to the 
velocity vector. 

This form (|5.2| ) shows clearly that the norm is not affected by that con- 
tribution to the equations of motion. Indeed, the partial derivative —777 

(i) 

ensures that the whole contribution vanishes when integrated over p). . We 
have the same structure that for the contribution (7.15) of [jjj or for the 
operator < 11 [0]|G| 11 [!(/)] >(°' 2 ). 



Since the first propagator in (5.2) and the second one in (^3) cannot be 
resonant, the ie can be dropped from them. 

We do not yet analyse the possible divergence of these contributions to 
< 11[0]|6|11[1(/)] >( 1 ' 2 ). 

5.3 Evolution of the distribution function 

We consider anew the case in which the two particles j and j' are perfectly 
localised with a well defined momentum (f4.7|) . If we perform the trivial 
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integrations, due to simplified form of the distribution function, we get: 
^ll[0](q (1 \p (1 W 2) ,P (2) ) 



ei 



1 1 



E^/a [11 /^ E E ^ 



(2vr)3 4vr^ 2 ., .. . a=1>2o=±1 



dp 



kW.v^-ibWy \ v ie + (il + akW).vO') - |l. V 0') - afcW 



x5(q^ - q J -)tf(q u ' J - - Pj )5(p 



H / dv\i 



or 



f v 0'). e («)[l]) e -*2-(%-<lj') 



(J) 



(5.4) 



^/ 11[0 l(q (1) ,P (1) ,q (2) ,P (2) ; 



(2») 



i=i,2 



a=l,2a=±l 



jfeW 



Mm 9 



ie + (Al + akWj.vO') - il.vCi') - ofcW / VakW.vOQ - ofcW 



xl 



Ov 



00 



,(a)[l] e -4.(q^q y ) 5(q (i) _ ^(qO") - q^ 



x£(p^ -pjO^P ' -P 



(5.5) 



5.4 Radiative reaction force 

The effect of the coupling of the Coulomb interaction with the field is to 
provide a supplementary force, the radiative reaction force, that changes 

the mean value of the momentum of one particle. The expression of the r 

(i) 

component iv of the radiative reaction force can be obtained by consid- 



ering the relation Fr = A < p\ >. By a partial derivative, we get that 

(i) 

the r component iy is provided by minus the coefficient of the expression 

is removed and where the variables 



„C?) 



( |5.4|) when the partial derivative 



dp 



T7T 
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qC?) ; qO' ) > p(i) anc i p(i ) are replaced by their values obtained from the Dirac 
delta functions. 



< >,- •• ! ' :: 



t (27r) 3 47r^' ej 



a=L,2 a=±l 



1 e -4-(«l3— Qj') 



kW.Vj-ifeW 



ie + (ii + akW). Vj - il.v,-, - ofcW j 



(PjuPjv) 
l~ 3~ 



(v,eWW), 



(5.6) 



(2 



a=l,2o=±l 



ie + (|1 + akW).v(i) - ±l.vO'') - a^W J VakW.v^ - afcM 



, (a)[l] 



dpi 



(5.7) 



We focus first on the first contribution. The value of that radiative reac- 
tion force depends on the relative orientation of the vectors position and 
momentum. We explicit the summations over u and v. We get then, using 
the value of the metric tensor g st to replace g st Vj s e[ a ^k^ by — v.e( Q )W k^ 
and g st v ]S e^ m k^ by -v.kNe£ a)[l1 : 



f\ 1 e^e 



< r >/_ '(2vr)3 4vr 
xpW-v.-.kWjeWN+Vj.eWI 1 ]^] 

(PjuPjv) 



a=l,2a=±l 



£ 



1 3 

(■mj+pj)2 (mj+pj)^ 

2 



(v,eWN) 



k [i]. Vj _ fc [i] y ^ e + (ii + a k[i]). Vj - |l. Vj v - aAiW 



(5.8) 
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The reality of this expression can be checked by considering the symmetry 
a — ► —a, 1 — > —1. 



5.5 Emitted power 

We distinguish the components of the radiative reaction force in the direction 
parallel and perpendicular to the velocity Vj of the j particle. The power 
emitted is given by < F^.Vj >. As can be seen, the magnetic force, arising 
from — g st Vj s (ef^ k*ft — e^' 1 ' A;] 1 ' ) does not contribute. The force parallel 
to qj — qji provides a radiative correction to the Coulomb force that is not 

considered here. We use Z^a=i,2( v j- e ^' 1 ')( v j- e ^ a ^ 1 ') 
obtain 



2 (v.ktW)^ 
V 3 (fe(W)2 t0 



< 



F^.v, 



>i= -t 



1 



e 3 e 3 



(2vr) 3 4vr 



-q,') 



a=±l 



(vj-k^W)- 
(fed 1 !) 2 

\ 2 / 



LkW ^ (l.p j )(p J -.kr 1 ]) 

(m 2 +p])2 



(m 2 +pj)2 



kW 



V; 



2^ v i' 



afcW 



(5.9) 



That expression is further analyzed in Appendix D, particularly in the sit- 
uation where the particle j' is much more heavy that the j particle. In the 
referentiel in which the heavy particle is at rest at the origin of coordinates, 
we treat the case where the vectors and Vj are orthogonal (the orbital 
situation). In such a case, all integrals can be performed explicitly and the 
final result is 



< F (j) .v,- > Iorb = -z&ev 



1 



o 3 3 , 2 2 . 3 3 /-. 



< F^.Vj >Horb= 7' 



rri; 



3 1 1 3 
1 1 



2\3 
3 I 



In 



1 



l + v. 



+ 



2vi 



(5.10) 



(5.11) 



9 3 3 ( 2 , 2\2 n 3 v 

These expressions enable to determine the component of the self electric 
field at the localisation of the particle: 



< E (Oj)- v j >orb= e)ej' 



1 
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-3 + Wvi 



3vj 



3(1 



,2^3 



2,v j 1 + Vj 



(5.12) 



For the geometry chosen, the radiative reaction force is known exactly by 
an explicit expression. 



5.6 Non-relativistic limit of the emitted power 

The previous expression can be developped in powers of u| to make the 
connection with the well known result. We have to consider the expression 
up to order u|. The result is: 

<E^(q J ).v i >^=|fe4 (5-13) 

3 rrij qj 

The coulombian acceleration of the charge j is provided by the dynamical 
function — £ = — J ,- v , . The mean value of its time derivative, due to the 
free motion ( |2.9[ ), is 

4 < *i >= / A» /<j 3 p« aszg™ ( -^>» > 

x£(qW-q^(p (j) - Pj ) 

In the geometry where the vectors and Vj are perpendicular, we have 
then ^ < a c j >= e ^' g 3 J • Therefore, we get the form 

< F^.Vj >7Vi ?0 r6= ^ej < — a cj > .Vj (5.15) 
If we restore the dimensions, we get 

< F^.Vj > NRIorb = — < a cj > - v i (5-16) 

The usual result, with a front factor |, is recovered directly, without having 
met any divergence for that contribution. This result is not astonishing. 
In the usual approach, the divergence appears as an infinite self-mass cor- 
rection, in front of the acceleration vector. Since we have considered the 



40 



geometry where the Coulomb force (hence the acceleration) is perpendicu- 
lar to the velocity, that divergence has no influence on the emitted power. 
In the general case, the contribution < F^) > n would provide the usual 
divergence. 



6 Emitted field due to the Coulomb interaction 

6.1 Subdynamics operator 

The determination of the emitted field due to the Coulomb interaction re- 
quires the determination of the creation operator in the first order in both 
the Coulomb and the transverse field interactions. The starting expression 
is: 

(1,1) 



< ll[l( ai )]|S(t)|H[l(/)] > K p'f 

— 1 r' 

= / dze~ lzt 

2ni J c 



z - kO').vO) - kO'). v (j') + fcW(mL 1] + m [ j) 

x (-i) e . e .,zL fd 3 l-l (— ^—)e\am~adl T )) 

{ 1)6363 2^ J d 1 P \dpU) dpU')) 

X {z- kO').vO') - k0''). v v ) + fcW(mL 1] + m™)) 

d 



dp? 



^ (v W.eC«W)(2-^-4 T (mgl-6) 



expbl k 1 ° ° 



1 



(6.1) 



z- kO').v(i) - kO').vO v ) + fcM(mj, iJ + m^ J ) 

for the order PF and a similar expression for the order FP. For the order 
PF, only the first propagator, at the extreme right, corresponds to a vacuum 
state while in the other order, FP, the first two propagators satisfy that 
condition and have to be considered inside the path c. We can proceed in a 
straightforward way, as in §5. 
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Since the operator A ([O]) can only deviate from unity when two field in- 
teractions take place, the expression of < ll[l(sj)]|S(0)|ll[l(/)] >( 1 ' 1 ) fn(f\] 
can be identified with the corresponding term < ll[l(sj)]|C'|ll[l(/)] >( 1,1 ) 
/rY(/)]- As the equivalence conditions imply /np,^.)] = /n[i( Sj )], the distri- 
bution function for the emitted field at first order in the field interaction 
and first order in the Coulomb interaction is determined in this way. 



6.2 Transverse emitted field 

We have therefore all the elements to deduce the emitted field (for sharp 
locations and momenta for the particles) 

<E^(x)> e ^ 1 ' 1 )= />kW E E fcW i e-(kW)ex P {m[kW.x]}e|e i ^ 

a=l,2 a=±l ~" 

Sli 1 1 ^ 



(2vr) 3 J I 2 \ie + (1 + akW).Vj - l.Vj-/ - afcM / (+ak[ 1 l.v j - ofcW 
+ / d 3 kW V V fc[ 1 ^e°(k[ 1 ])exp{m[kW.x]}(-) e 2 e7 v^i 



a=l,2o=±l 
(i 3 / 1 



(2^) 3 7 I 2 \ie + (1 + akW). Vi - l.vj/ - aJfeM J \akW.Vj - ofcW / 

J.) 1 2 ^|^ c (a)[l] e -i(|l +O kW].« b -l.q J 0. (6 . 2) 



This new expression is the equivalent of ( [4.8D in presence of the Coulomb 
interaction. It determines the field due to the accelerated particles in terms 
of the actual values of the position and momentums pj of the charged 
particles. Usually, expressions of the acceleration fields are given in terms 
of the retardated positions. We look for the comparison only for the radia- 
tive force, since we have illustrated in §4 the equivalence of the formalisms 
outside the locations of the particles. The self- field of the particle due to 
the Coulomb interaction, is then given at first order by 

< EHaj) > e ^= -eh'lj^s J rf3k[1] E E <( k[1] ) 



a=l,2a=±l 
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I 2 We + (1 + okWj.Vj- - l.vjv - ofeW J (+akW.v j - afcM 



, l v 



2 1 1 

+6jej, ^(27r)3 



X / (f 3 /- 



(vj.e 



| d 3 kW ^ E <( k[1] : 



a=l,2a=±l 
1 



I 2 \ie + (1 + akW).Vj - Lvy - afcM / ^akW.Vj - afcM 

, ^£ e (a)[l] e -il.[q,-q 3 .,] 



(6.3) 



Using (3.14), that expression can be identified with the result obtained from 
the operator (in the previous section) that leads to the usual expression 
for the self-force in the low velocity limit. 



7 Finite classical electrodynamics. 

7.1 General outline of the approach 

We intend to prove that resummations, usual in the context of statisti- 
cal physics, enables to get rid of the divergences in the kinetic operator, 
computed from the expression of the self-field at the location of the particle. 
These resummations involve a renormalisation of the propagators associated 
with the particles. Our analysis proceeds through several steps. 

In the first one, the previous divergent contribution for the self-field is 
written in an adequate way. A linear dependence in a cut-off wave number 
K c is found, as expected. A usual diagrammatic representation, such that it 
can be found in Q for instance, is convenient. The dependence of a simple 
cycle on a cut-off wave number K c is then analysed and established. The two 
propagators of the previously divergent contribution are then renormalised 
using simple cycles. When their dependence in K c is reported in the ex- 
pression of the self-field, it induces a supplementary and the resulting 

c 

contribution vanishes! 

In the second step, we use cycles that are renormalised by themselves: 
the propagator present in the expression of the cycles is a renormalised one. 
That first resummation contains only disjoint or inserted cycles in arbitrary 
numbers, with no overlap: two vertices of a cycle are either disjoint or 
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completely inserted with respect to two vertices of another one. The con- 
tribution of the cycle is thus defined in a self-consistent way. The resulting 

I 

dependence of the renormalised cycle in K c is non analytic, in K 2 ■ This re- 
sult of the self consistency can be understood in a simple way. The integral 

defining the bare cycle behaves as K c . The renormalised cycle involves one 

_ i 

propagator that, by hypothesis, provides a factor K c 2 . When combined, 

i 

their product reproduces correctly the assumed K 2 dependence. When that 

value is introduced into the expression of the self-field, each of the two prop- 

_ i 

agators provide a factor K c 2 that compensates the K c dependence due to 
the integration. Therefore, the limit K c — > oo is well defined. 

That property of the renormalised cycles can be extended to all contri- 
butions: For imbricated vertices, each new wave number will be associated 
with two propagators, that can be renormalised by (renormalised) cycles. 
The dependence of these new vertices with K c is then the same as the the 
simple cycle. Therefore, they can also be used for a global renormalisation 
while preserving the finiteness of the self-field. 



7.2 The previous divergent contribution 

The correlated elements describing the self-field are easily obtained from 
the expression of the creation operator. We have to focus on the divergent 
contribution to the the self-field. Its computation starts from the expresson 
of < ll[l( Si )]|S(i)|ll[l(/)] >( 1 ' 1 ) (fltUD) and we have to reconsider it in 
order to be able to add the contributions of higher order that will ensure its 
finiteness. 

The final dependence in the cut-off wave number K c has to be appre- 
ciated after the a last vertex, enabling to get or the self-field from the 
creation operator, has been introduced. A displacement operator involving 
the wave number exchanged through the Coulomb interaction is natural and 
will remain in our expressions. We keep the dependence in particle operator 
that has provided, in the previous "II" contribution, the ultraviolet diver- 
gence, in the computation of the self-interaction. We still restrict ourselves 
to the effect of the absence of fields: all strenghts of the fields oscillators are 
close to zero: the distribution function of the action variables of the field 
oscillator are peaked around zero (the limit is taken afterwards). 

We reconsider the contribution, affected by a d-index, responsible for the 



divergence, from the first steps of its computation. In the expression (|6.l|) 
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the operator 1. ^ Q ® U) — 7^77) has to act on — ^(v^.e^M). Moreover, in 

the factor f 2— m-( m fl' ~~ &) I > onr y the derivative will provide a non 

vanishing contribution for the self-interaction. 



< 



ll[l( Sj -)]|S(t)|ll[l(/)] >i u) 



= -/' 

2m Jc 



dz e 



-izt 



k(i).vO') - k0 v ). v 0') + ^^(mL 11 + m' lh 



Maifi) ifcU 7 )) 



l 



z - kO').v(i) - k(i'). v 0'') + fcW(mL 1] + mj 1 ), 
1 



x H) e .?7TV E E 



1 ._^ y , (v ,, e » m) ) 2 _^ 



/? J 



expfoi -kW. 



5 



9 



9k0) dm f 



kO).vO') - kO''). v (i') + fcWtmjJ 1 + mJJ, 1 ) 



(7.1) 



That expression has to act on the distribution function /[Y(/)j describing 
the absence of field. We change the place of the displacement operators (in 
variables) to obtain: 



< ll[lfo)]|S(f)|ll[l(/)] fj {f)] 



x— f 

2-KZ Jc 



, n N 3 E E 

(2vr)2 p =lj2 b=±i 



i \ 2 



exp?J -kW. 



8 



9 



9k0) 9m W 



e 



-izi 



(kC?) + bkWj.vC?) - kO'). v (j') + bkW 



x( i)ejej 



-1 

2^2 



3 



TrCv^.e^W)') 



d 



dr]f _ 



z - (k(i) + 1 + 6kl 1 ]).v(i) - (kO'') - l).v(i') + bkK 
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- (ktt) +l).vW - (kO v ) -i). v O'i, 
xeKafr-sMffoN - ^(^ - ^O^/^M o- (7-2) 

The previous creation operator, enabling the computation of the self-field 
is obtained from that expression by picking the residue at the pole corre- 
sponding to the last propagator. 

<^(*) ^M^^ / ^ E E <(kW) 



(27 

1 / 1 



Q =l,2a=±l 

x I d 3 l" 



'I 2 \ie + (1 + akW).Vj - l.v^ - ofcW / 



akW.Vo — afeW / dp™ 



(7.3) 



Each propagator provides a decreasing factor while the jacobian to get 
to spherical coordinates provide a (fcM) 2 factor, hence the linear divergence. 



7.3 The simple cycle 

We intend to renormalise to propagators through simple cycles in a first step. 
In order to be acquainted with their contribution, we explicit the operator 
< ll[0]|£(z)|ll[l(/)] > (0 ' 2) fX, f)] that contains only the cycle. The two 
interactions involve the same particle j. 

We consider only contributions that will not contain displacement op- 
erators involving the wave number of the interacting field mode acting on 
the distribution function. Indeed, in their presence, a convergence factor for 
ultraviolet modes would arise from the dependence of the involved distribu- 
tion function in that wave number. The selected terms act as c-number in 
the particles variables, hence the c- index. Those contributions are the one 
responsible for the divergences since the other contributions involve char- 
acteristics of the distribution function. In the contribution ( |3.15| ), we have 

thus b = —a and — ttt has to act on (v^.e^W). We act on the field vacuum 

dpi 3 ' _ y ' 

distribution function /[Y(/)i to get: 

< n[o]|£(*)|ii[i(/)] >?' 2) /w/)] 
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(2tt): 



[1] 

'2 



E ^,o^,o E E & [[^^-^(^^ 



a=l,2o=±l 



dp. 



(j) 



exp a 



-kW. 



a 



[1] 



^ _ k(j). v Cj) - k0''). v v ) + fc[i](mL 1] + m [ lY 



(-i)e 



7 3 

J (2tt) 5 



E uiT 



/3=1,2 



_ I(V (,»,(«,( 2 ^ __a (m | t l + 



exp(-a) { -kW. 



5k0) 9m W 



z - kO'). v O') - kO'O.vO'O + fcl^mj 1 + mjj, 1 ^ 



Jl] 



(7.4) 



When acting on the field in its ground state, to obtain a non vanishing 
contribution, we need an interaction with the same mode (a = (5). The 
summation-integration over field variables can then be performed: 



< ll[0]|TC(z)|ll[l(/)] >f s > 

{ Z7T ) J a=l,2a=±l 

1 



z-k0'). v 0) -kO'O.vC?') 



x 



[fc[l] e (-)[l] _^ v 0) (e HW4l] _ e (a)[l] fell])]] 



kO').v(i) + akW.vC?) - kCiO.vC*"') + fcW(-o) 



-2tt 



a 



^ v 0). e H[i]) 



-kW.vO') -kO'O.vO'O), 



(7.5) 
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We define the cycle contribution C(z,v^) by 



^ J a[l]=l,2a=±l ^ 



1 



+ akW.vW) + fcW(-o) 



-2vr 



9 



U) 



0) 



(vW.e 



(7.6) 



and we recognise C(z — k^.v^ — k^ ' ).v^ '\ v^) in ( |7. 

Let us examine the dependence of C in the cut off wave number K c . 
The first matrix element contains a contribution finite for A;M — > oo. If 
we take into account the summation over a, the propagator will involve a 
(&M)~ 2 dependence. The whole contribution behaves therefore linearly as 
K c . If z approaches the real axis fom above, as it is required from the 
construction of the subdynamics operator, we can use z = y + ie and upon 
replacement in ( |7.6| ), the expression of C(i/ + ie,v0)) contains a Dirac delta 
function, of argument y + ok^.v^ + A;M(— a) and a principal part. The 
contribution due to the Dirac delta function is independent of the cut off 
while the contribution due to the principal part provides the K c dependence. 

Moreover, it is easily recognised that C(z,v^) vanishes for z — ► 0. 
Indeed, in that limit, the contribution that could arise from the Dirac delta 
function <5(akM .v( J ) + £;M (—a)) vanishes since the vanishing of its argument 
cannot be satisfied (except for fcM =0 but the jacobian provides a factor 
and the contribution due to the principal part vanishes by parity in 
the summation over a. 



7.4 Renormalisation by simple cycles 

We first consider the consequences of inserting a simple cycle in the contri- 
bution (17. ll) . For the sake of illustration, we consider the insertion of such 



a cycle before the vertex corresponding to the Coulomb interaction in (7.1). 
Details of computation can be found in appendix E. The insertion of the 
cycle involves of course a supplementary propagator (see the square on the 
first propagator in ( |E.3D ) and the previouly considered cycle contribution 
(7.6) with the same argument as the new propagator. That property holds 



for the insertion of an arbitrary number of cycles, enabling a renormalisation 
for the propagators. The resummation of the class of diagrams where all 
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the cycles are at the left of the particle vertex (Coulomb interaction) is: 

< ii[i( Si )]|s(*)|n[i(/)] >£ ,oo) fj (f)] 



H^tt'tt E E ( !I kY expb \~ k 

(2vr)2 / 3=i,2 6=±i V k 1 I 



d 



<9k0) dm^ 

— /' ciz e^ 2 * f z - (k^ + 6k). - k^.v^ + &£; 

27TZ Jc V 

C(z - (k^ + 6k) .v^ - k^").v^") + bk, v^)) _1 



-1 /• o 1 



1 



Z 2 V z - (kO') + 1 + 6k). vO') - (kO'') - l).vO'') + 6fc 

1 



(; 



(k(?) +l). v 0') - (kC?' # ) -l).vO'') 



[i] 



(7.7) 



When an arbitrary number of cycles are also added in the two other 
propagators, the contribution is: 



< ii[i( ai )]|s(t)|n[i(/)] >i 1,oo) fj m] 

1 ^ 



(-i)e 



J (2tt)2 ptt,2bt^l ^ k 



E £ f exp6 



x - J: /"' <fe e^ 2 ' f z - (k^ + 6k). - k^.v^ + 6fc 
-C(z - (k^ + 6k).v^ - k^.v^ + bk, v(^)) _1 



(z - (ktf> + 1 + 6k) .v^ - (k&"> - l).v^") + bk 
-C{z - (k^ + 1 + 6k).v^ - (k^ - l).v^ + bk, 



-i 



x (z - (k^ + l).v^ - (k^ - l).v^ 

-C(z - (k^ + l).v^ - (k^ - l).vk"\ v^)) _1 
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xe L (^-5kky s{v W _ m)5{r] M _ W )8^ ] q 5^ } q . (7.8) 



For the evaluation of superoperator, we have to compute the residues 
to the (multiple) bare poles due to the last propagator. By resummation, 
the procedure corresponds (cf. Lee model (2j) to the computation of the 
dressed poles. In view of our previous analysis, C(ie,v^^) = 0. A solution 
6 of 6 = C(0, v^) is thus 9 = 0. 

The solution of z- (k^ +l).v^ - (k&") - l).vk") - C{z- (k^ + l).v&') - 
(k^") - l).v( J ''),vW) = 0. can thus be written as z = (k^') + l).vW + 
(k^ ' — l).v^ the bare and dressed poles coincide. The residue at the pole 
disappears when the creation operator C is computed from the product CA. 
Therefore, for our purpose, we do not have to bother about a renormalisation 
of the vacuum propagator and we can continue to use the simple propagator. 

In place of the second term in (6.5) , we get the same expression with 
the cycle contribution in the propagators 

< E^(q,) >J^= e)e r 1 --^ f £ £ e?(kW) 

x J d 3 l^ (ie + (1 + akW). Vj - l.v jt - ak^ 
-C(ie+ (l + akW). Vi -l. Vj v -afcW, Vi )) _1 

J ! \ 

^akW.Vj - afcW - C(ie + okW.Vj - afcW, Vj) J 

x ^^£ e («)[l] e -il-[q,-q,-']. (7.9) 
upj v 

In the referentiels where the velocity vy vanishes, we get 

< E^-) >^ 1} = ^^77^3 / d3kW E E <(kW) / cft± 

" I 'v J a=l,2o=±l J 



1 



ie + (1 + okM). Vj - ofeW - C(»e + (1 + akW).Vj - afcM, v,-) 

: ! ) 

akW.Vj - afcM - C(ie + akN.v,- - aJfcW, v^) / 



x l 7 >)[l] e -l-[cb-q f ']. 



'/■',. ■ • (7.10) 

mj 
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In the low velocity limit Vj — > 0, the previous expression vanishes by sym- 
metry due to the angular integration, (the integral over 1 and k become 
independent but the final result could be divergent.) 

We introduce a cut-off function in ( [7.10 ). The same cut-off will be used 
in all future expressions. 

" \"") J a=l,2a=±l J 1 



ie + (1 + akW). Vj - afcM - C(»e + (1 + akMj.v,- - afcM, v,-) 



akW.Vj - akW - C(ie + akN.v,- - aJfcW, Vj) 



1 / r fa)ni r -»l.fq-q / l 



(7.11) 



where in the low velocity, neglecting the magnetic component of the force 
and simplifying the computation of the derivative with respect to pj , 



C(z,Vj) = -—ej 



d 3 k 



1 



z + k.Vj — k 



k 2 + K 2 ' 

We examine the behaviour for small velocity Vj. 

1 



C(z,0) 



16vr 2 „ 1 
-el 



rrij 3 (27r) 3 Jo 



dk k l 



z — k 



+ 



1 



z — k.Vj + k 

(7.12) 



z + k 



Kl 



k 2 + K 2 
(7.13) 



For Qz > 0, and $tz = y, we prove in appendix E that 

y 



C(y + ie,0) 



^ 2 1 r7ri ,3 



(2vr)S 



C (y 2 + K 2 ) 



+ i 



?Tij y 2 + if 2 



(7.14) 



We insert the value of C{y + ie, 0) in the expression (|7.1l[) , neglecting 
the dependence of C on the velocity Vj. That dependence is not seen as 
capital to get the qualitative behaviour while the remaining dependence in 
the denominator is required for the coupling between the integrations over 
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the wave numbers I and k. We have: 



7T (2TT) ., , 
v ' a=l,2 a=il 

/.ft 1 ' 1 U 1 



Z 2 ^te + (1 + okWj.Vj - ofeW / - afcM j 

[(H-okW).Vj- - afcW] 2 + *T 2 
X o 

[(1 + akW). Vj - afcW] 2 + K 2 - ^-e 2 ^ c 3 - i^^ 2 [(l + akM). Vj - afcM] 

X 1 ^ 

[akW.Vj - afcW] 2 + K c 2 - ^-e)Kl - i^K 2 [akW.Vj - ofcW] 

xi^.HLh,^,,.^. (7 . 15) 

We compare that integral ( [7.1 5| ) with the contribution arising from (6.5): 
< E^) >?™ = efe,,i * / ^ E E e?(kW) 



(27T) 

i / 1 

d 3 l 



a=l,2a=±l 

1 / 1 \ / 1 



I 2 \ie + (1 + akW).Vj - afcM / ^akW.Vj - afcW / 

(7 - 16) 

we note the presence of two similar supplementary factors. Except for the 

3 

domain of values where the wave numbers are of the order of K£ , the de- 
nominator in these factors are dominated by — —e 2 K 3 In the relevant 
range of values for the wave numbers, thanks to the cut-off function, the 
numerator in the factors are of the order of K 2 . Therefore, these factors 
can be replaced by to evaluate the dependence of ( 7.15|) with respect 



j 

to K c . Since the resulting integral is known to diverge linearly with K c , the 
resulting dependence vanishes as K~ . The introduction of the (simple) cy- 
cle produces a convergence with respect to K c that is much more too strong: 
we go from a divergence in K c to a convergence in K~ . The consideration 
of renormalised cycles will produce the correct behaviour. 
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7.5 The renormalised cycle C 



The renormalised cycle C is introduced by replacing, in the definition of the 
cycle, the free propagator by a fully renormalised one. C is thus defined by 
a self-consistent equation: 



c(zM j) ) 



(27T) 



a[l]=l,2a=±l 

' 1 



+ okW.vW) + 



-2vr 



9 



U) 



(J) 



(v^.e 



C(z + akW .vC») + jfeW (-a), 

if 2 
fcW 2 + K? ' 



(7.17) 



We analyse the behaviour of C(z, v^) with respect to K c and v^). The self- 
consistent equation ( |7.10| ) is diagonal with respect to the velocity v^) . In a 
first investigation, we consider the situation in which v^) is close to zero. In 



these circonstances, 



dp, 



™(v®.e^W) ) is — and the summation over 



r can be performed and we use the unity of the e^^ 1 ! vector. We get: 



C(z,0) = (-2vr)-e 2 ' 



x 



E 



m j (2tt) 3 
1 



m 2 + k? 



a=±l \ 

We have (C(z) = C{z,0)) 



C(z + A;[ 1 l(-a),0) / 



(7.18) 



C{z) 



(-47T) 
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(2tt) e 



fcW 2 + ^ 2 



1 



a Si U + ^ [1] (-a) - C(z + fcM(-a)), 

X 2 



16vr 2 2 1 P „ , 2 

f I ^ 

- fc - C(* - A;) 



+ 



m 2 + e: 2 
/ 1 



+ fc - C(z + fe) 



(7-19) 
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In our units in which the velocity of light is unity, the dimension of C{z) is 
the inverse of a time (as z) or the inverse of a wave number as k. We can 
define K m by ttt. We can therefore write C(z) = zc(-J—, S^) and we have. 



, z K c 

zc{ iT'ir j 



l 

1 



dk k z 



Ki 



k 2 + K 2 



'VH K,„ ' RV, 



+ 



z + k-{z + k)c(^^ 



(7.20) 



We call 7 the ratio between the cut-off wave number K c and K m . Introduc- 
and u = jr—, we have the equation for y real (from above): 

,2 



c(y,i) 







duvf 
1 



7 



7 2 + u 2 



+ 



\(ie + y-u)[l- c(y - u, 7)] / \(ie + y + u)[l- c(y + u, 7)] , 

(7.21) 

What kind of consistency can we deduce from ( [7.2l[) about the dependence 
of c(y, 7) with respect to 7? 

Let us be more specific about the behaviour of c(y,j). We assume that 



c(y,7) = «7 r + (5 



(7.22) 



in the range \y\ >> 7 and the self consistency should determine the complex 
parameters a, (3 and the real parameters r and s. The parameter r is 
assumed positive and s has no definite sign a priori. That expression for 
c(y,7) asssumes that lim^oo c{-y\) = ^" ^ e ^ rs ^ cnec k that property in 
appendix F by computing the difference A(y, 7) = c(y, 7) — c(— y, 7) for large 
positive y. The asymptotic behaviour of c(y, 7) in the domain y >> 7 >> 1 
is then checked in a self consistent way by the proposal ( 7.22| ). In order to 
make predictions on the behaviour of the self-field when the contributions 
of the renormalised cycles are included, we need mainly the behaviour of 
c(y,7) in the other domains, in particularly 7 >> y >> 1. 
From the expression 



c(y,7) 



1 

y Jo 



+ 



du 
1 



7 



7 2 + u 2 



1 



(ie + y-u)[l-c(y-u,j)} 



(ie + y + u)[l- c(y + u,-y)] 



(7.23) 
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and our previous results, we have to deduce first the behaviour of c(y,j) in 
the other domains: y « 7 and y in the range of 7. 
We assume that 

c(y,7)=7^(-) (7.24) 
7 

in the range \y\ « 7 and the intermediary range. Self consistency should 
hold and determine the complex function <?(-). g(0) is assumed 7^ while 
g(^) is assumed to vanish for ^ — > 00. That expression is analysed in details 
in appendix F and self consistency established. 

That dependence can be further used to evaluate the dependence of the 
electric field with respect to the cut-off value K c (7 = j^-). 

7.6 Behaviour of the self-field 

We can now look at the consequences for the self- field. The computation of 
the residue to get the required creation operator can be pursued in a similar 
way. The residue at the vacuum pole k^.v^) + k( J ).v« ' does not play a 
role for the creation operator: only the value of the pole matter. We recover 
therefore an expression similar to the expression ( |7.11| ) in which the simple 
cycles are replaced by the renormalised one's: 

< <M >7& A) = ^77^3 / E E <(kW) 



a=l,2o=±l 



I 2 \ie + (1 + akWj.Vj- - ofeW - C(ie + (1 + okWj.Vj - afcM, Vj 
1 



akW.v,- - ofcW - C(ie + akW.v, 




The same manipulations as for the simple cycle lead to a global dependence 

_i 

as the product of a factor X c arising from the integral and two factors K c 2 
due to the presence in the denominators of the renormalised cycle function 
C. Therefore, the previous expression ( 7.25| ) for the self electric field pro- 



vides a finite result in the limit K c — * 00. We note that the dependence 
of the self-field on the charge ej disappears also in that limit. This prop- 
erty is not at all astonishing and prove that an effective cut-off for wave 

2 

numbers larger than \ is present. The final dependence is proportional 
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to that value and the dependence in ej cancels out. Therefore, we meet 
here a typical case where a non-analycity in the expansion in the charge in- 
duces divergence on individual diagrammatic contributions but an adequate 
resummation enables to obtain finite results. 

7.7 Beyond the inserted cycles 

If we consider two cycles with an overlap, with respect to the previous 
case, we get two supplementary propagators. If the propagators are renor- 
malised by the previous considered inserted cycles, the contribution of the 
imbricated cycles implies a possible 7 factor arising from the new integra- 

tion over a wave number and two Ki factors arising from the denominators. 
Therefore, the global behaviour of the imbricated cycles is the same as a sin- 
gle (renormalised) cycle. The imbricated cycles can thus be introduced for 
renormalising the propagators and the previous analysis is still valid. There- 
fore, if the contributions are computed in terms of renormalised propagators 
with all possible vertices (defined also in terms of renormalised propagators), 
the result of the resummation is a contribution to the self-field that is finite 
in the limit of an infinite cut-off K c . The replacement of a bare vertex in 
the expression of the self-field by dressed vertices does not change either the 
finiteness of the global contribution. 

8 Conclusions 

Our present work have illustrated the feasability of a reformulation of clas- 
sical electrodynamics, that takes explicitly into account the corrections due 
to the self-fields. Moreover, the procedure avoids the existence of runaway 
solutions: causality is an ingredient of the construction of the subdynam- 
ics operator. Therfore, our expression for the self-force is not in terms of 
the time derivative of the acceleration but involves the actual position and 
velocity of the charged particle. We justify in that way the procedure pro- 
posed by several author to avoid the runaway solutions: the replacement of 
the time derivative of the acceleration by the time derivative of the external 
force. In the traditionnal approach, the self-force is naturally computed from 
the characteristics of the trajectory and the replacement has to be added 
by hand. Here, we have made the opposite step: our expression in terms of 
the mean field has been shown to be equivalent with the usual expression in 
terms of the time derivative of the acceleration. 



56 



The present approach constitutes a statistical description of interacting 
charged particles and electromagnetic fields: we are far from the classical 
view of well defined values for the variables associated to the fields and the 
particles: all these variables are statistical with a joint distribution function 
that evolves with time. The use of a reduced formalism enables to treat 
the distribution functions that are the most relevant for the computation of 
mean values of all the dynamical functions. 

Two distinct ingredients are required. The first one is a relativistic sta- 
tistical description of interacting fields and charged particles in which no un- 
observable potential appears as dynamical variables. Balescu-Poulain have 



developed further the ideas of Bialynicki-Birula [ 21 1 , |2^] and his coworkers 
to provide such a formalism free from dynamical constraints. The elimina- 
tion of the Lorentz condition is a key element of the present work that avoids 
the usual derivation of the self-forces via the Lienard-Wiechert potentials. 
The second ingredient is the possibility, that we have developed in collabora- 
tion with C. George, of getting rid of the self-field by defining an appropriate 
subdynamics. When both elements are combined, we obtain a finite kinetics 
for the description of the interacting charges and fields in which no explicit 
self-energy process is allowed: the kinetic operator takes into account all the 
effects and its computation, althought lenghty, is straightforward. 

The formalism developed in this paper offers the basis to tackle in a 
new way the divergences in classical electrodynamics, through (infinite) re- 
summations of diagrams. The renormalised propagators are then defined in 
terms of them-selves (for a similar procedure, see |24|) and the solutions of 
the resulting non-linear equations do not admit a simple expansion in the 
charge, enabling the presence of the natural cut-off, linked to the classical 
radius of the electron. 

The present paper illustrates only one of the multiple potentialities of 
the approach. Many problems can be aborded within the present formal- 
ism, such as the charge renormalisation, for instance, of higher order effects. 
Moreover, we have considered the charged particles outside an external in- 
fluence: the distribution function corresponding to the field vacuum has 
been used thoroughly in this paper. The effect of the magnetic field has not 
been specifically considered: when computing the power dissipated in the 
motion, its effect disappears. We have not taken advantage of the statistical 
nature of the formalism: a sharp distribution function has been assumed for 
the positions and velocities of the particles. A statistical nature for the field 
has also been ignored. 

An irreversible extension of CED, analogous to the treatment of the Lee 
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model in quantal case, requires the construction of the generators of the Lie 
associated with the extended dynamics. The relevance of such an extension 
is still to be established. 

We dare to state that we have presented new tools for dealing with 
problems in classical electrodynamics. The approach may look tedious but is 
nevetheless straightforward. It opens new perspectives, not only in classical 
electrodynamics but also in quantum electrodynamics [^(|. Obviously, a lot 
of work remains to be done. 



9 Appendix A. The Lorentz force 



In this part, we consider only one particle interacting with a free transverse 
wave. The particle will be pointlike, with a specific value for the velocity. 

From the expression (|3.13|) for < 11[0]|6|11[1(/)] >( 0,1 ), we deduce the 
change to the one particle distribution function due to that contribution, 
assuming the independence of the field and particle variables. 



dtf(Kp,t) 



1 



(2tt)5 



d 3 k^ / dri[ 



4 1] 



E 



W a=l,2a=±l 



E E 



x exp a 



)[i] 



9 st v, 



. \yt K r e r H )\ q 



-kW. 







a 



^ dm, 



m \ f(k,p,t)f m ( n f\m [ ^rj [ ^m [ ^M 1] )- 

(A.l) 



If we suppose that / describes a particle localized at some place r(i), 
/(k, v, t) is proportional to exp— ik.r(i) (|3.6|) . The action of the displace- 

'' ^ can thus be performed and we get 



ment operator exp a 
easily: 



dtf(Kv,t) 



m 



(2tt)! 



d 3 k^ / dri[ 



4 11 



E E 

a=l,2o=±l 

{»kW.r(t)} 
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x /(k, v, t)f w (r?! 11 , -ad a>1 , r§ , -a5 a , 2 ; k™ ) . (A.2) 

The mean values < E^(x) > and < B^(x) > of the fields can be deduced 
easily from (U) and (HJ): 

<e- ( x) >= / d 3 k [i] r dv m r dv m f\in f 1 ^ 

J Jo Jo Jo Jo 

x- 2 ^ E E ^^(k^r/JCkW) 

(27TJ2 Q=l 2a= ±l 

x exp{m[kW.x - 2vrd 1] (k [11 )]}/li](x [1] ; kW) 

=/>kw /°v /v 

J JO JO JO JO 

x^^ E fc [11 ^e r Q (k [11 )r ? l(k [11 )exp{ia[k[ 1 ].x-2^ a (kW)]} 

(27r)2 a =l,2a=±l 

[i] [i] 

1 /* /"OO /"OO 1 1 

= _L_ / d 3 k [l] / / ^1 V V fc[^e°(kW)r,I(kN) 

(2vr)W Jo Jo Q=li2fl=±1 

x exp{mk [11 .x}/[ 1 ](r/f ] ,-a^ jl ,4 1] ,-a^a,2;k [1] ), (A. 3) 



< b^x) >= / d 3 kW / 4 11 / 4 11 / ^S 1] / de. 

J JO JO JO JO 

x-^TT E E fc wi (-l) a e?(km)uJ(kW) 

(27TJ2 a= i i2a =±l 

x exp{m[kW.x - 2n^K^ ] )}} f[l](x {1] ^ 

d w r drj w r dr)? f di m 

Jo Jo Jo Jo 



[1] 



_L_ £ £ ^ 1 ^(-l) Q e r "(k[ 1 ])r ? l(k[ 1 ])exp{m[k[ 1 ].x-2< Q (kW)]} 

(27r)2 a= i i2 a=±l 

^M^)^, , mf 1 , 4 1] , mf 1 ; k™ ) 



[i] [i] 



1 /* pOQ /"OO 1 1 

(2tt)2J Jo Jo £1,2 £±1 

exp{mkW .x}f {1] (77W , -a5 a ,i, j£ ] , -<w5«, 2 ; k [1] ) , (A.4) 
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so that we can proceed to the identification (331). 



10 Appendix B 

This appendix completes the list of the matrix elements of < 11[0]|G|11[1(/)] >(°' 2 ). 
< 11[0]|G|11[1(/)] ><M) 
2 1 d 



x [ju[l] p (a)[l] _ ^,(i) (>)[%[!] _ f WWf.W 



\orja' T]W , 

E5 [i] n 5 [i] _ exp 2a < — . — —. .. 



1 



-akW.vC?) +afcW 
d 



d 



dm, 



[i] 



+ £(- 



i a 



(2tt)3 ^0) 
fc [i] e (a)[i] - 5 ^')( e f)N4i] - e W[%fl) 



or 



■akW.vO') + afcM 

E 5 m w/ m w,o ex P 2a |- kW - ° l> 



[1] [1] 



1 '^kO') ^ 



(B.l) 



< 11[0]|G|11[1(/)] >gf 



E fW 1 f 4 1] E E f?IV (IT 

s % o£) J Jo 11 Jo a % a hA m ) \ m 



i 



fc [i] e (a)[i] _^ v p-) (e («)[i] fc W -eWWfcf) 
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E ^Aw.o^P - 2 « k[1] - 



3 d d 
a — — a- 



fikW 5mL 1] 0*& ] 



+ E 



i=i,2 



_L d_ 

(2vr) 3 apO") 



1 



-okW.vC?) + afcW 



W\ 5 /Jl]\ 5 



[fc [l] e ( ? ')[l] _ /t '^) (e (a')[l] fc [l] _ ^IjJlpJ) _ 2afc W 
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[i] [i] 

,m 2 



" a 



(B.2) 



< 11[0]|6|11[1(/)] >g) 
1 3 

X 



5 / [IK \ 

\ 1%, I 



akW.vO') -afcW 



E ^ [1] [1] n eX P i _a ~~7T 



(?) 

r' 



a' 



+ E(- 

J=l,2 



fc [i] e H[i]_^0) (e WW4i]_ e H[i] fc W) 



-z e„- 



(2 
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n [i] mi '° ™ 2 '° I SmL 1 dm% J 



[i] 

,m 2 



11 Appendix C 

The expression of the complete electric field < E(x) > e j(°' 1 ) ( |4. 12 ) is eval- 
uated explicitly, using its identification with < E _L (x) >a j< '°' 1 ''- Multiplying 
numerator and denominator by (fcM _|_ kM.Vj), we have 

< E 1 (x) >?' (<M) = (2*) ej pipj | <i 3 klH Binlkl'lfx - q,)] 

* ( (*M)»-Uv,)- ) ^ +k ' 1 '^' ) " fc ' I,V ') W 
Let us place the x axis along (x — ) and the y axis along v±j , defined by 

v,.-v [v J .(x-q J )](x~q j ) 
V ^ - V J ^qJF • 

eCnil 1 f+°o r+oo r+oo 

<E (x) > a ] ' = (47r)ej / / cZA^ / dA^ sin k x |x - q 3 1 

J —oo J — oo */ — oo 



1 



A; 2 - ( 



x ((k x e x + k y e y + k z e z ) - kvj) -. (C.2) 

The contribution involving /c 2 e 2 vanishes obviously for parity reasons. The 
integrand has to be even for a simultaneously change of the sign of k x and 
k y . Therefore, 



e (0 1) 1 f+°° r+oo r+oo 

<E (x) > a 3 ' =ej— ^ / dk x / dfey / d/c 2 sin£; x .|x - q,-| 

J —oo J —oo J —oo 



e 



J r+oo r+oo r+oo 

''I 
1 -) 

+ Vjyky) J 

1 

1 r+oo r+oo r+oo 

j~o~2 / ^ / ^ / dk z sin k x \xi- qj\ 

<27T J —oo J —oo J —oo 

{{k x G x + kyGy) — {vj x k x + v ] jyky)~v j) . (C.3) 



yfc 2 {Vj X k x + Vjyky)^ J 

x (k(k x e x + kyBy) — k(vj x k x + 



A^ (Vj X k x + Vjyky)^ 
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We use dimensionless variables of integration. We then replace sink x by 
±r(e lkx — e~ lkx ) and perform the integration over k x by residue at the pole 
of -rj—, — ttt — tyi in the correct half plane. We have 

1 \ 



r+oo r+oo r+oo 

h = I dk x dk v / dk z e ikx 



-oo J — oo J — oo 



{Vjxk x + Vjyky)^ J 

x [(k x e x + k y e y ) - (v jx k x + v jy k y )vj] . (C.4) 
The pole is obtained by the equation 

k {Vj X k X + Vjyky) = 0, 

V jx) ~ 2k x kyVj x Vjy , 

Therefore, 



- *4) - ^xkyVj x Vjy + - ) + ^ = 0. (C.5) 



for — 



k y v jx v jy ± ^ (kyv jx v jy ) 2 - (1 - w] x )[A;2(l - u^) + / 

kyV jx Vjy ± 1^(1 - Vfj^l - + fef] - (kyVjxVjy) 2 



(C.6) 



Due to the factor e* fca: , the relevant pole for ii corresponds to the plus sign 
and we have 



/+oo r+oo 
dk y I dk z e 
-oo J —oo 

1 



^J^-^ x )[kl{i-V 2 jy ) + kl]-(kyV jx V jy f . kyV jx V jy 



kyVj x Vjy + «y/ (1 - - ^) + _ {kyVj x Vjyf 



(1 " t£) 



hOO /" + 00 



y/(l-^ x )[kl(l-^. y ) + kl\-{kyV jx V jv ^ .kyVj^y 

= ir I dky I dk z e (1 v i x> e 

— oo J —oo 



:i-v] X Wy(l-v] y )+kl]-(kyV jxVjy y< 
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kyv jx vjy + iJ (1 - «y[^(l - v) ) + A; 2 ] - (k y v jx v jy ) 



+& y [e y - u^vj]] . (C.7) 
We replace the oscillating factor according to its parity in k y . 

1 



ha = ni 



(1 " v%) 



2 \ l e x V 



/+oo r+oo 
dk y / dk z 
-oo J —CO 



xe 



cos 



KyVjxVjy 



(C.8) 



Introducing polar coordinates r and # in the k y , k z plane, we get 

r2w 



/*oo rz 

I la = wi[e x - v jx \j](l - v 2 jx ) drr 



d6 



rVj x Vj y cos I 



— r./l— v 2 —v 2 cos 2 9 i 
Xe V jx jy cos ( 

X r°° r 

= -iri[e x - v jx Vj](l - v) x ) j drr J 



d6 



-r[y/l-v'j x -vj y cos 2 9+ivj x v jy cos 0] _|_ e ~ r l\/ /l ~ v %~ v ' A j y cos2 9-iv jx v jy cos 6] 



(C.9) 



The integration over r is readily performed. 

r2w 



1 /" 2 

ha = ^i[ e x ~ V jx \j\{l - V 2 jx ) 



dB 



^ l - v2 jx- v2 jy cos29 + iv ^ v 3V cos °] 
1 

[\Jl-v] x -v] y COS 2 ~ iv ixV jy COS #] 2 ^ 



= me. 



r 

-v jx -Vj](l- v] x ) 

(i 



2?r 1 - f %. - v\. COS 2 

d(9 3 - m 



V ]x v jy cos2> 



I 1 - *4 - v2 iv cos2 + cos2 e ? 



1 /" 

= iri\e x -v jx Vj\ T -j- / 

- v jx ) Jo 



jx "jy 

'? 2 — 7? 2 
jx u jy 



2w 1 — V 2 rr . — V?, COS 2 ' 

d6- 



jx"jy 
v 2 x v] y cos 2 1 



[1-^cos 2 ^] 2 



(CIO) 
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Taking cp = 29 as new integration variable, we get 

1 



'la 



n a a. 1 ~ V i* ~ Hl/C 1 + U P ~ hjy( l + v j J cos ( 



[1 



(C.ll) 



From formulae 2.554.2 and 2.554.2, 148 of p9|, we read 



A + B cos x 



1 



+ 



(aB — ^46) sin x 

(a + bcosx) n (n — l)(a 2 — b 2 ) [(a + 6cosx) n_1 
( Aa - 6J3) (n - 1) + (n - 2) (a£ - A6) cos x 



(a + b cos x 



in— 1 



-dx 



A + B cos x 



B Ab-aB 

-dx = T x + 7 

a + b cos x b b 

with, formula 2.553.3, for a 2 > b 2 
1 . 2 



1 



dx, 



a + b cos x 



-dx 



b 2 



: arctan • 



a + b cos x 
V 'a 2 — b 2 tan § 



(C.12) 
(C.13) 

(C.14) 
(C.15) 

(i-t&)(i-t&)S 

We now turn to the second term of ( p.7|) that is evaluated in a similar way: 



a + b 



Therefore, the last integration can be performed and we get 

„2l 



2 2 [1 — IHj 

/la = 27r - v jx )e x - Vj x Vj y e y ] 



1 



'16 



2?r 2 i- 



1 



3 v jx v jy ^ _ ^2 



[1 ~ «?)] 



(C.16) 



"3V> 3X1 

The sum of the contributions I\ = l\ a + lib is the contribution along e x of 
Ii a ( p.!5| ) and is given by 

h = 2n 2 z[l - v]}- 1 -— T e x . (C.17) 

The contribution from I2 is obviously its complex conjugate and, from 
and ( C.4 ), we have for < E _L (x) >a the expression: 

1 



< E (x) >a 



[1 - v 2 ] 



ej(0.1)_ «, 1 1 o 



J 2vr 2 2i |x- qi | 2 
1 1 



/1 



(C.18) 
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12 Appendix D 



We evaluate first in this section the power dissipated by the radiative force 
< F^.Vj >i Q5.9| ). The second contribution is treated afterwards. 

We decompose the vector kW into its component k , and perpendicular 

k[j^ to the velocity vector Vj. The scalar product (l.kW) becomes the sum 
(l.lcj-| +l.kJJ ). By symmetry, the last term will generate a vanishing contri- 
bution when integrated over kjjL The remaining scalar product (l.kj, 1 ^) can 
be written as pJ 2 (l.p_j.)(kW.p 3 -) and combined with the other contribution. 
Since -K — , T , = $, J. we get 



< v, >z= / J i £ 



(27T) 
l. Pi )(p i .k[ 1 



a=±l 



fcWy Ue + (il + ak[ 1 l).v i 




(D.l) 



1 e^, 



(2vr)3 4vr p 2 (m 2 +p 2)§ 



d 3 k [1] i <in jl 



xe 2 



%.(q J -q/)( Lp .)(p.. k [l]) 



•PjAPj- 

2 



kW.Vj-JfeW 



^ + (il + kW).v i -il.v i /-A;[ 1 ] 



^ie + Cil-kW^-^Lv^ + ifeW 



(D.2) 



We can consider a situation where the particle f is much more heavy that 
the j particle. In the referentiel in which the heavy particle is at rest at the 
origin of coordinates, we have: 



1 e]ey 



in- 



(2vr)3 4vr p 2 (m 2 +p2) § 



<fk ] - I d 3 l± 
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;e - i J*(l.pj)(p J .k['l) 

.( ! 

\ie + (li-km). Vj + km 



ie + (ii + kW).v J - m ) 

(D.3) 



We consider first the case where the vectors and Vj are orthogonal (the 
orbital situation) . We place the x axis along q., and the y axis along Vj . We 
have: 



1 e$e r 



//•CO POO POO j 
d 3 /^/ / / cfl^cfl, 
J -co J — oo J -co l x ' L y ' l z 

xe- il -fV{l yPj ){ Pj lS) 



v ie + (i/j / + 4 11 )^ -^W 



ie + (ii v -4 1I K + *; [11 



(D.4) 



The integration over Z y can be performed by residue, closing the path in the 
upper plane Ql y > 0. Indeed, the integrand decreases at least as The 
only pole to be considered is l y = + 



<FW.v J>w = 1 ^ 



m; 



J 



d 3 k^ r r 

J — CO J- 



x(i^l 2 x + l 2 zPj )( Pj kW) 



(2tt) 3 4tt pj(mj+pj)% 
dl dl 2m e'^V 



^e + (^ v ^T7f + 4 11 K-^ [11 . 

1 

ie + (^yfT7!-4 1] H + ^ [1 



(D.5) 



The ie can now be dropped since they have played their role in determining 
the relative position of the poles in the complex plane. We introduce polar 
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coordinates in the l x , L plane to get 



< F^.vj > Iorb -- 
2m 



1 e]ey 



(2tt) 3 4vr p 2 (m 2 +p 2)| 



1 



dll 



2tt 



d0 



- k$)vj + km 



Since 



(D.6) 



*(ky_uj — k — ) — ^ we have 



1 e| ej v 



m 2 Apjf 



(2tt)3 4 p 2 (m 2 +p2) f 



2- 



xe" 



j|gjCose^,[l] 



2(4 1] ^-/« [1] ) 



d0 



(D.7) 



By definition, Jq W d6 cos (y cos 9) = 27rJo(y), do being the Bessel function. 
Therefore, 



< F^.Vj >i 



orb~ 



4tt ejef 



(2vr) 3 4 p| (m 2 +p2) 



d^W 



dll 



(i^.)2 + ^\ v . _ k \x]y 
From p. 686 of [19] we have (formula 6.565.4): 



(D.8) 



J v {bx)x 



dx 



o (x 2 + a 2 Y+ l 2^r(/x + 1) 



(D.9) 



where K u (z) is a bessel function of imaginary argument (— 1 < 'Rv < K(2^ + 
|), a > 0, b > 0). We can apply that formula for x = I, with i/ = 0, /x = 0, 

b = ^Qj, or = -^- a — ^ — . The function i^oO 2 ) is represented in 8.432.1 by 



the integral {y = 0): -K"o( z ) = fo° e zc dt. The integral over I can thus 
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be performed: 

4 . 1 if0( # 1 -^) ) , (D , ) 

We take fc? 1 = cos 0, x = cos 0, J d 3 kK . . . = J™ dfcW (fcW) 2 J+ 1 J 2 *" d<\. 

<F - Vj>/orfc - (2vr)3 4 ^( TO?+ ^)H^ } ^ 
x-^— K ( g ^ M(1 -^ ) ). (D.ll) 

Xfj — 1 f j 

The formula 6.561.16 p. 684 of @ is: 

r ^(ax)cfa = 2^ 1 a-^ 1 r(i±4±^)r(i±^), (D.12) 
JO 2 2 

with + 1 ± i>) > 0, 5Jia > 0. That formula ( p,12| ) can be applied for 
fcW, with = 2, i/ = 0, a = Mlz^il. 



:>, \ 2 e^ej/ m 2 



7rV* V 2V 4 ( m 2 +J ^)|t,2 

x / +1 — f ^ ) 3 , (D.13) 
J -i xv j — 1 \qj(l — xv j) J 

J A^V 4 ( m 2+^)l^y-i (1-x^) 4 

(D.14) 

The last integral is direct and leads to: 

4 ?7i^ i?^ 1 

< FW. Vi >/ orb = -efc, J ^ (D.15) 

In the other geometries, some integrals are not known explicitly but can 
be shown to be more convergent that the orbital case that provides a finite 
result. 
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We now turn to the second contribution. In place of ( |D.4| ), we have now 
(by a change of variables, 1 in this expression is of ^1 in the I contribution): 



< F (qj).Vj >Horb= e j e j'~ 

1 



1 1 



K + (p0))2)l 



Pj 



oo roo poo 



Pj v j 



7T (2tt) 3 

(4 [1] )fa-4 [11 

(fed 1 !) 2 
{ Pj k y [l] ){v 3 k^ ) 



(fed 1 ]) 2 



OO J — oo •/ — oo 

1 



dl x dly dl z 2 



(D.16) 



Vie + (Zj, + 4 1] )vj - « + (iy - 4Vj + / 

The integration over Zy can be performed by residue, closing the path in the 
upper plane Ql y > 0. Indeed, the integrand decreases at least as Z~ 3 . The 
only pole to be considered is l y = iy/l 2 + I 2 . Introduce polar coordinates in 
the l x , l y plane, we get 



1 



Vjin 2 



<F- (qj ).v J >^=e3 e ,_^ + ^ ))2)| 



d^ 1 



X 1 



(A; 



V ) 



1 



2 7! 



d6 e~ ilxq i 



{il + ky 1] ) Vj -k^J V (il - 4 1J ) v j + k[1] 



,4\-&w 



(D.17) 



Since 



fying the Jo Bessel function in d9 cos (ycosO) = 27rJo(y), we get: 

,2 



-2(4%--fc [1] ) 

(%) 2 +(4 1 'fi-fe [1] ) 2 



and identi- 



< F (qj)- v i >iiorb= -e,-e 3 -/ 



1 



x 1 



(4 [1] ) 2 ' 

(fc(W) 2 



dZ Uo(lqj) 



(2tt) 2 ( m 2 + (p(j))2)| 

2 



(z^-) 2 + (4% - 



(D.18) 
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Using (p.S|) and formula 6.561.16 p. 684 of IS] leads then to: 



< F ± (q j ).v j > IIorb 



7T 



r£) 



2 2 
V 3 m 3 



1 



( m 2 + (p0-))2)l q] 
\ 3 



f+1 / I 

' I dx (1 - x 2 ) - - 

-1 K^-XVj) 



(D.19) 



The last integral can be performed to provide the result ( |5.11| ) of the main 
text. 



13 Appendix E. The simple cycle 
13.1 Insertion of a cycle 

The vector kM in ( [7,ip . is replaced by k with a corresponding change in 
the other variables associated with the field mode. To avoid ambiguity, rjs 
are replaced by ep that has the same use of deferring a vanishing limit. We 
have to act on the field vacuum to use the previous result: 

<ll[l( Sj )]|E(t)|ll[l(/)] >% 3) fJ if)] 

= zl f dze -^ ( I 

2vri Jc \z- k(j).vO') - k&'O.vW) + k(rn a + m a >) 



afll=l,2 o=±l 



[jfe[l] e («)[l] - g«vV\4 a)[1] kW - e^k?)) 



kW.vO') + akW.vO') - kO'O.vO'O + jfeW(-o) + A:(m a + ^ 



-2vr 



3 



(v^.e 



x 



1 



z - kO').vO') - kO'O.vO'') + k(m a + m a /) 



i 



kO').vO') - k(j').v<>'') + k(vn a + m a / 



J (2vr)l 



E E 

/3=1,2 6=±1 
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exp 6 < — k. 



3 



8 



dkO) dm, 



z - kO'). v O') - kO'O.vO') + ^(m, + mp>) ^ 

■Kr srKr 



Usual manipulations lead to: 



< ll[l( S ,)]|S(t)|H[l(/)] >S' 3) fj {f)] 



(-i)ej 



-1 



(2tt)! 



E E 

0=1,2 6=±1 



VP 



a 



—izt 



/ / eXp6 \- k -5kO) 9m, 



1 



X 2iriJ c dZe \z- (kC?) + 6k) . vO') - k(i') . vO'') + 6A; 



1 



(k(j) + 6k).v(») + &£: + okW.vC?) - kO'O.vO') + fcW(-a) 



— 1 /" 1 



z - (k(j) + 1 + 6k). vO') - (k(i') - l).v(i') + bk 



(E.l) 



From the definition of the cycle, we have 

1 



z - (kO') + l).vO') - (kO'O - l).vO v ) 

(E.2) 



< ii[i( ai )]|s(t)|n[i(/)] >S' 3) fj 



(2tt)2 



0=1,2 6=±1 



A; 



9 



X 27riJ c dze \z- (kO) + 6k) . vtf) - kO') . vO") + 6fc 



-izt 
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xC(z - (k^ + 6k) .v^ - k( J ").v( J ") + bkM j) ) 

-l r .* . 1 / i 



(k(i) + 1 + fck).v(i) - (k(i') - l).v(i') + 6fc 







l.^ By ,(v«.e(»)j 2- 

1 



(kO') + l).v(?) - (kO'O - l). v O'0 



xe 



(E.3) 



13.2 Evaluation of the simple cycle 

For Qz > 0, and 5ftz = y, we have 



%C(y + ie,0) 



16vr 3 



1 



(2tt) 3 7 



(ffcfc 2 - k) + 5(y + k)\ 



2e 2 K 2 y 2 



rrij y 2 + K 2 



k 2 + £T C 2 
(E.4) 



Let us note that the sign of 9C(y + ie, 0) is unusually positive and QC(y + 
ie, 0) vanishes at y = 0, in accordance with the solution 2 = of z— C(z) = 0. 

For very large y (y >> K c ), it behaves as a constant 3 m : c . Since the 
imaginary part of C(z) vanishes for z = 0, in the denominators in ( |7.9|) and 
( |7.10| ), that property does not involve a new ie rule in the computation of 
the subdynamics operator. For the real part of C(y + ie,0), involving the 
principal part, we have 



3?C(y + ie,0) 

16vr 3 , 1 
-ej 



rrij 3 (27r) 3 jo 



dk k' 



K 2 

c 



k 2 + K 2 



(y-k) 



(y + k)' 



167r3 2 1 , T r ^ 
-e?T^r(/l +/2) 



m. 



(2vr) 5 



(E.5) 



= lim 

R^oo Jo 

lim 

R— >oo 



l 



l 



ln(R -y) + ln(|y|) + ln(i? + y) - ln(|y|)] = 0, (E.6) 
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dkKi 



M I dkK 2 



1 



(y-k)' y (y + k)- 

1 , 1 

+ 



K 2 C 



k 2 + K 2 
K 2 



(ie + y - k) (ie + y + k) 
1 Ki 



k 2 + K 2 



roo 

' -oc dk Tie -n,- k)k- \ l\- 



ttK 



II 



C (y 2 + K 2 Y 

I2 vanishes linearly for y —* and behaves as 7ri£f | for large y. 



(E.7) 



14 Appendix F. Self-consistency analysis 

14.1 Consistency check on A(y, r y) 
From the expression ( 7.21| ), we have: 



A(y,7) = c(y,7) - c{-y,j) 



y Jo 



du u 



2 r 



7 2 + u 2 



1 



+ 



1 



, (ie + y - u) [1 - c(y - u, 7)] / \(ie + y + u)[l- c(y + u, 7)] 

!- 1 t. — : — *w 



(ie-y - u) [1 - c(-y - u, 7)] / \(ie - y + u)[l - c(-y + u, 7)] 

(F.l) 

As for (2.30), we separates the contributions arising from the 5 functions 
from the contribution arsing from the principal parts. That separation does 
no longer correspond into a separation between a real and an imaginary 
part. The contribution Ag(y, 7) from the 5 functions is: 

, 2 



1 f°° 7 

y Jo 
1 



7 2 + u 2 



x U(y-u) 
+8(y + u) 



+ 



1 



[l-c(y-u,i)] [1 - c(-y + u, 7)] 
1 1 



+ 



[1 - c(y + u, 7)] [1 - c(-y - u, 7)] 



1. . 1 2 Y 

y [l-c(0,7)] 7 2 + r 



(F.2) 



As(y,j) contributes to the second term of ( 7.22] ) since it decreases as -. 



The parameter s is thus related to the behaviour in 7 of 



[l-c(o, 7 )J • 
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The contribution A-p(y,j) from the principal parts is: 



Ap(y,7) 



x \V 



1 



1 r j 2 7 2 

- / duu 

y Jo 7 + 

1 



(y - u) 
l 



(y + u) 



[1 - c(y- it, 7)] [1 - c(-y + u,7)] 
1 1 



[l-c(y + u,7)] [l-c(-y-u,7)] 



(F.3) 



To determine the behaviour for y >> 7, it is too simple to say that we have 
two y factors in the denominators in that expression and thus a behaviour 
as 4j-. We have indeed to make sure that the convergence of the integral 
is not affected by the limit. Obviously, if we neglect u in front of y in 
that expression, the remaining integral over u diverges linearly. We analyse 
therefore the behaviour of the integrand I(u). For the first contribution, 
the form ( 7.22|) cannot be used since the argument y — u of c(y — it, 7) is 
not large for all values of u inside the domain of integration. For that term 
Api(y, 7), we separate the domains < u < 2y and u > 2y in the conditions 
y >> 7 » 1. In the first domain, the domain of integration over u is finite 
and we have: 



Apia (y, 7) 
xV 1 



U* 2 T 2 
duu 



y Jo 



7 2 + u 2 



y-u 



[1 - c(y - u, 7)] [1- c(-y + u,j)] 



(F.4) 



1 c^y ( 1 

Apia(y,7) = — / du-y 2 V 

y Jo \y- 



1 



1 



[1 - c(y - u, 7)] [1 - c(-y + u, 7)] 

1 f 2y , 7 4 ^ / 1 
H — / an^ rrT 7 



y Jo 7 2 + n 2 \y — u 



1 



.[1 - e(y- n,7)] [1 - c(-y + u, 7)] 
The first term A-pi ao (y,7) can be written as: 



(F.5) 



7 f y /I 

Apiaa(y,7) = — / dvV[- 

y J-y \v 



[l-c(-v,j)\ [l-c(v,i)] 



(F.6) 
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where the principal part symbol can be dropped. For large y, we can replace 
±y as limit of integration by ±00, providing the convergence of the integral. 
For large v, 



1 



L[l 



V 



[1 - c(u,7)] 



1 



1 



[v(l - Q7 r ) + PY] 



[v(l — aj r ) — P"Y S 



[v 2 (l-aY? + P 2 l 2s X 



(F.7) 



That factor behaves as i and the integrand in ( |F.6j ) behaves as \, ensuring 
the convergence. A-pi aa (y,7) contributes also to the second term of ( 7.22| ). 

The second term A-pi a t>(y,Y of ( |F.5| ) behaves as least as 4j since no 
convergence problem can arise. 

For the second contribution Ap2{y, 7) of ( |F.3| ), the asymptotic behaviour 
(|7.22|) can be used in all the integration domain and we have, keeping the 
two terms: 

,,2 



Ap 2 (y,7) 



1 

y Jo 



du 



7 



7 2 + u 2 



V 



y + u 



1 



Elementary manipulations lead to: 



Ap 2 (y,7) 



l 

y Jo 



duu 



T 



-2PY 



Y + u 2 [(y + u) 2 (l — cry 
^VikHil) can be split as (cf. ( |F.5|) ): 

Ap 2 (y,7) = —[ du 



r~\2 



2~,2sl 



/? 2 7 



(F.8) 



(F.9) 



7 

y Jo 



dm 



y Jo 
1 



u) z (l — 07 
-2/?7 s 



r\2 



2„,2sl 



/? 2 7 : 



7 2 + u 2 [{y + n) 2 (l — 07 



T^2 



2-,2sl 



/? 2 7 



(F.10) 



In the first contribution, the integral converges and behaves as ^, providing a 
global behaviour at least as In the second term, the factor -z+^p provides 
an effective cut to the values of u that contributes in the integral. We can 
in the integrand neglect u with respect to y and the global behaviour is in 
p-. Therefore, the form ( |7.22j ) is compatible with our analysis of A(y,7). 
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14.2 First consistency check on c(y,j) 



We now turn to the analysis of c(y, 7) itself. We can choose a positive sign to 
y to fix the ideas since A(y,7) can then provide the behaviour for negative 
y. We start with (cf. (BO)) 



c(y,i) 



1 

y Jo 



00 'V 2 

2 7 



duu' 



7 2 + u 2 



(ie + y-u)[l-c(y-u,j)] 



1 



(ie + y + u)[l-c(y + n,7)] 
For y > 0, we have 



(F.ll) 



cs(y,i) = — {-ni) 
y 



1 



cfot u £(y — u) 
7 2 



7 



y y ^[l-c(0, 7 )] y 7 2 + y 2 ' 
The contribution c-p(y,j) from the principal parts is: 



[1 - c(y - it, 7)] 7 2 + n 2 

(F.12) 



+v 



1 



du u z 
1 



7 2 + v? 



V 



y-uj [1 - c(y - u,7)] 



y + u/ [1 - c(y + it, 7)] 



(F.13) 



For the first contribution to c-p(y, 7), the form ( [7,22 ) cannot be used since the 
argument y— u of c(y—u, 7) is not large for all values of u inside the domain of 
integration. That integral cannot be written as a sum of integrals involving 
different principal terms since only the global integrand in (F.13) decreases 
enough to ensure the convergence at infinity. The following decomposition 
still holds (10) : 



7 2 f c 

y Jo 



+v 



du 
1 



V 



1 



1 



y-uj [1 - c(y - u,j)] 



y + uj [1 - c(y + u, 7)] 



+— / du 

y Jo 



+V 



7 2 + u 2 



V 



1 



1 



y-uj [1 - c(y - it, 7)] 



y + uj [1 - c(y + u, 7)] 



(F.14) 
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For the second term, the factor 



ensures the convergence and cut the 



integral for values of u limited by 7. For that part, the dominant behaviour 



in y is obtained by neglecting u in front of y. For the first term in (F.14), 
the convergence in u is ensured by the other factors. For studying the 
convergence, we take u much larger that y and we can replace the c by their 
value ( 7.22| ) and the integrand I\(u,y) is 



h{u,y) = 

y 



1 



1 



y — u [1 - ay 



0^- 



h(u,y) 



.11 

y 



+ 



+ 



1 



1 



y + u [1 - aY 



0-% , 

^ y+ui J 



(F.15) 



(y — u)(l — «7 r ) — 0j s (y + u)(l — aj r ) — 0y 



h(u,y) = ^2[y(l-aY)-07 s ]^ 7T 
y u z [l 



(F.16) 
(F.17) 



a i' r ~ [2/(1 — «7 r ) — 0y 

Therefore, the factor that ensures the convergence of the integrand in u con- 
tains as a factor y. The integral for op(y, 7) contains therefore a contribution 
independent of y for large y. It corresponds to (and determines) the first 
term of ( |7.22 ). The 7 dependence factor in the y independent contribution 
for very large y is thus 7 2_r , that should be identified with 7**. The self 
consistency possibility is thus r = 1. 



14.3 Second consistency check on c(y, 7) 

In the expression ( 7.23|) for c(y,7), in the right hand side, three domains 
for the u integration variables can be distinguished. In the first one, the 
argument of c(y — u, 7) and c(y + u, 7) are inside the conditions asssumed 
<< 7. In another domain, we are in the conditions studied in 
\y\ >> 7. We have also the transition domain where 



for ( |7T24D : \y\ 
the preceding section: 
y is of the order of 7. 

We separate anew in (7.24) the cs(y,j) and cp(y, 7) contributions. If we 
consider y > 0, we have still the form ( p?.12j) : 



cs(y,j) = -(^rr 



7 



oc 



iriy 



c(0, 7 )] 7 2 + y 2 c(0, 7 )" 



(F.18) 



Assuming the form (7.24) in cs(y,j) provides a behaviour as yj 2 as dom- 



inant contribution, much smaller that the assumed ( 7.24 ) behaviour 
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For the second term c-p(y,j), we can use ( F.13| ): 
1 f°° 7 2 f / 1 \ 

y Jo Y + u 2 l Vy- «/ 



y- uj [1 - c(y - 11,7)] 

+V {VT^) [l-c( y + n, 7 )]}- (F ' 19) 

Combining the behaviour ( |7.24| ) and the behaviour ( 7,22j ) where r = 1, we 
assume the following behaviour for ep(y, 7) for large 7. 

op(y )7 ) =7^(-)+75i(-), (F.20) 

7 7 

2 

where the function 51 has the qualitative features of ^jW^p is a very 
large number): it is negligible for y « Ij and becomes 1 in the other limit 
y » I'y. This expression reproduces qualitatively the behaviour of ep(y, 7) 
for all values of y with respect to 7. To check the consistency, we introduce 
that expression in (4.4) for y « 7. We have (1 is negligible in front of 7) 

^ 1 r , 2 7 2 \„( 1 \ 1 



7 y ^0 7 + u \ \y-uJ M 




( —J—") — * }> . (F.21) 



If we introduce variables t = ^, v = ^, we get 



t - v 



[7 2 y(* - v) +J9i(t - v)] 



+V ( — 1 (F.22) 



The front factor 72 can be simplified and we have 

9 ^~~lIo dUU 1 + v 2 { V ( T^v) [g(t-v)+ 75i (t - vj] 

+V (i + W [g(t + v) + 7gi (t + v)]} ' (F ' 23) 

The remaining dependence of 7 plays a role only for very large values of t — v 
or t + v. Its consequence is to diminish the contribution for that part of the 
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domain of integration. Since that expression has no problem of convergence 
due to the factor 1 * a , that dependence does not play a sensitive role in 
the evaluation of g. It can therefore qualitatively be dropped and we can 
deduce that c-p(y,j) behaves as 72 for the the relevant domain of values of 
its argument. 
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